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A Computer Code for a One-Dimensional 
Dynamic Model of the Mesosphere and 
Lower Thermosphere 


1. INTRODUCTION 


T.J. Keneshea, one of the authors, pioneered studies of the numerical solution 
of sets of continuity equations in the atmosphere. He was the first to find that these 
equations were ''stiff'' and were not solvable efficiently by any numerical technique 
then available. His efforts!»2,3 led toa computer simulation of the composition of 
the upper atmosphere that could be run inexpensively on the large scale computers 
available in those days. After these computer codes were applied to the undis- 


4 


turbed atmosphere~ and to several natural? and manmade atmospheric disturbances, 


it became apparent that they could not predict long term atmospheric behavior, 
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1. Keneshea, T.J. (1962) A Computer Program for Solving the Reaction Rate 
Equations in the E Ionospheric Region, AF CRL-02-035. 
2. Keneshea, T.J, (1963) A Solution to the Reaction Rate Equations in the Atmos- 
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phere Below 150 Kilometers, AFCRL-63- 
3. Keneshea, T.J. (1967) A Technique for Solving the General Reaction Rate 
Equations in the Atmosphere, AFCRL-67-0221, AD654010. 


4, Keneshea, T.J., Narcisi, R.S., and Swider, W., Jr. (1970) Diurnal model of 
the E-region, J, Geophys, Res, 75:845-854. 
5. Keneshea, T.J., and Fowler, R.J. (1966) Computed Electron, Ion, and 


Neutral eos Profiles for the Solar Eclipse of 12 November 1966, 


Such simple computer codes, in which only chemical reactions among the atmos- 


pheric gases are considered, are excellent for determining short term variations 
of the atmosphere, but give completely erroneous results for sclutions beyond a 
few hours. In fact, these codes are valid only for solution times shorter than the 
transport time constants for the various atmospheric gases. 

It then became obvious that any realistic atmospheric simulation should be 
based upon the equations of motion for the particles as well as their continuity 
equations. This report describes in detail a computer scheme for solving large 
sets of coupled partial differer.:ial equations of conservation of mass and momen- 
tum and also presents the computer code itself with instructions for using it. Also 
inciuded are the results of an application of the code to the development of a one- 
dimensional midlatitude winter simulation of the composition of the mesosphere and 


lower thermosphere. 


2. THE EQUATIONS OF STATE 


Tre equation for the conservation of particles of a particular gas within a mul- 


ticomponent gas mixture is 


On; ; ag 

me Pah eve, (1) 
where n, is the concentration of species i, P. is the rate at which it is formed, and 
nL; is the rate at which it is removed through gas phase chemical reactions with 
other members of the ensemble, The vector flux of the species, 9;; is induced by 
a combination of all nonchemical processes, In this report, only molecular and 


eddy diffusion will be considered in the flux term 


$. =n.C, + n.Vv. (2) 


where GC: is the velocity of species i resulting from molecular diffusion and Vv, is 
the velocity resulting from eddy diffusion. Since only the vertical distribution of 


species is to be considered here, the continuity equation can be written as 


at = P, - nL; - ae (aw yt n.v;) (3) 


where w, and v, are the vertical components of the velocities con and vi respectively, 


The equation of motion for a neutral gas particle in a multicomponent gas mix- 


ture can be written as& 


oii gig gs - (Cc. -G,) (4) 
ar ae ae 


where Pp; * nm, is the mass density of species i, a; is its thermal diffusion factor 
in a multicomponent mixture, P; is its partial pressure, T is the temperature of 
the mixture, fis the acceleration of gravity, and Vij is the collision frequency of 
Species i with species j, The equation of motion is written here in terms of the 
collision frequencies between particles rather than the molecular diffusion coeffi- 
cients in a multicomponent gas as was done by Keneshea and Zimmerman’ and 
Keneshea et al.8 The pressure of the gas is simply 


p, =n, kT (5) 


where kis Boltzmann's constant, Therefore, in the one-dimensional case 


1 P; kT | 1 On, (1+@)OT 1 
Dy .+a@,.—VinT + g = — |-- + ——— e+ 


where m, is the mass of species i and Hy = kT /mg is its scale height. The equa- 
tion of motion in the vertical direction for a neutral gas particle can then be written 


as 


Ow, kT 1 On, (1 +a,) OT 1 
at om weet eat a Dow (7) 


If inertial effects are ignored, then the equation of motion for ions and electrons in 


the presence of both magnetic and electric fields are? 


6, Chapman, S., and Cowling, T.G, (1970) The Mathematical Theory of Non- 
Uniform Gases, Cambridge University Press 

7, Keneshea, T.J., and Zimmerman, S,P, (1970) The effect of mixing upon 
atomic and molccular oxygen in the 70-170km region of the atmosphere, 
J. Atmos. Sci, 27:831-840., 


8. Keneshea, T.J., Zimmerman, S,P,, and Philbrick, C.R. (1979) A dynamic 


model of the mesosphere and lower thermosphere, Planet. Space Sci. 
27:385-401, 


“, Banks, P.M., and Kockarts, G, (1973) Aeronomy, Part A and Part B, 
Academic Press, New York, a 


x +—Vp. -F- = —C, =-lv,.(C, -C, 
at 7,7 e a (E ee C; x B) ; Y,; ( i C;) (8) 
and 
ec, 1 1 a ae 
+ ~p-—(E+—C B) = - Cc -C 
Sf oe te ae ee 2s V6 | 27 ©) 


where q is the particle charge, m is its mass, E is the total electric field vector, 
c is the speed of light, and B is the magnetic field vector. After much algebra, 
ignoring the horizontal component of the electric field as well as the horizontal 
velocities of the particles, and assuming that the electron and the ion temperatures 
are the same as the temperature of the neutral gas mixture, the equation of motion 


for the vertical velocities of the ionic species becomes 


Ow. kT 1 On. 1 On 2 OT 1 
a ae ee £ 
ot m, {ny OZ Ne OZ T OZ Hy 


2,/z 2 
Pea a Sore Lar 
=-w. Vp a (3 ‘) + Yo viw; (9) 
i j jog 
J Q? ce we 2 j 
1 J 


where n, is the concentration of the electrons, I is the magnetic dip angle, and 


(10) 


where By is the magnitude of the magnetic field and e is the charge on the electron, 
The eddy diffusion velocity, taken from the formulation of Colegrove et al, 10is 


1 On, 10T 1 
wi KS a2 tr ag tT 


where K is the vertical eddy diffusion coefficient and H is the scale height of the 
mean mass, 


10, Colegrove, F,D., Johnson, F.S., and Hanson, W.B, (1966) Atmospheric 
composition in the lower thermosphere, J. Geophys. Res, 71:2227-2236, 


10 
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Whenever a photodissociation rate coefficient is needed in the calculations, the O5 
and Ox column densities are estimated and a simple table look-up is made for the 
particular a required, 

A short discussion of each photodissociation process and the photodissociation 


ratu coefficient tables are presented in Appendix B. 


6.2 The Photoionization Rate Coefficients 


Reactions 185 through 191 and reactions 206 and 207 in the list of reactions in 
Appendix A are the photoionization processes considered in this report. The EUV | 
solar fluxes, the total absorption cross sections for Op, No, and O, and the ioni- | 
zation cross sections for Oo, No, O, H, and He are listed in Table B21, The EUV 
solar fluxcs from 1 to 121.57nm are from the compilation of Manson, oe The X-ray | 


fluxes below 1nm are taken from Swider., Ae The cross sections for total absorption 


and for total photoionization of Oo: Nos and O were obtained from a consideration 
of reported experimental and theoretical values. References, definitions, anda 
discussion of the absorption spectra can be found in Huffman? for wavelengths to 
30,.4nm. For shorter wavelengths, the values given by Swider?? were used, A 
discussion of some photoionization rate coefficients can be found in Keneshea and 
Huffman, 21 The photoionization cross sections for N, H, and He are taken from 
Banks and Kockarts, 2 

The photoionization of the metastable 0,(7A g) molecule has been included in 
these calculations. The wavelength region of interest is from 102.7 to 111,.8nm, 
where ground state O» is the primary absorber of the incident solar radiation but 
is not ionized itself. The effective absorption coefficient for Oo and the ionization 
cross section for O,('Ag) are taken from Huffman et al.?? Nitric oxide is ionized 
readily by Lyman o@ radiation, The ionization cross section for NO at Lymana is 


2x 19718 cm” as given by Watanabe. ae 


18, Manson, J.E. (1976) Satellite Measurements of Solar UV During 1974, 


AFCRL-TR-76-0006, AD A021490, 
lo, Swider, W., Jr. (1969) Ionization rates due to the attenuation of 1-100A non- 
flar:; solar X rays in the terrestrial atmosphere, Rev, Geophys, 7:573-594, 

20, Huffinan, R.E. (1972) Photochemical processes. Cross-section data, Chap- 

ter 12, in Defense Nuclear Agency Reaction Rate Handbook, DNA 1948H. 

21, Keneshea, T.J., and Huffman, R.E. (1972) Solar Photoionization Rate Con- 

stants and Ultraviolet Intensities, AFCRL-72-0667, 

22. Huffman, R.E., Paulsen, D.E., Larrabee, J.C., and Cairns, R.,B, (1971) 
Deoercase in D~region 09 (1A4g) photoionization rates resulting from CO, 
absorption, J. Geophys. Res. 76:1028-1038, 

Watanabe, K. (1954) Photoionization and total absorption cross sections of 


vases, I. Ionization potentials of several molecules, Cross sections of 
NH 3 and NO, J, Chem, Phys. (22:1564-1570, 


te 
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The total rate coefficient is obtained simply by summing Ji(A) over all the wave- 
lernsths, 


-% F(A Ni 


= 1 
J. 5 (AV TA) e 


J (54) 


6.1 The Photodissociation Rate Coefficients 


Reactions 192 through 205 in the list of reactions in Appendix A are the photo- 
dissociation processes considered in this report. In computing the transmittance, 
only molecular oxygen and ozone are considered as absorbers of the solar radia - 
tion, In addition to these two species, the absorption and subsequent photodissoci- 
ation of H,9, H,Oo, NOp, NOg, N,O, HNO), COp, and HO, are discussed in Ap- 
pendix B, Also presented in Appendix B are the photodissociation of N,O;, HNOg, 
CHy, and CHO, even though these species are not included in the aeronomic cal- 
culations discussed in this report. 

The incident solar flux TA) and the absorption cross sections for Oo, 03; 
H,9, H5Oo, NO», NOg, and N,O used in calculating the photodissociation rate 
coefficients are listed in Table B1. The cross sections for Oz in the wavelength 
region of the Schumann Runge bands (175-205nm) are not tabulated here. In this 


wavelength region, the molecular oxygen dissociation data of Hudson and Mahte?> 


are used, The incident solar fluxes listed here are those tabulated by Ackerman, ae 
In the wavelength region 123 to 190nm, however, the Ackerman fluxes have been 
replaced by the measurements of Heroux and Swirbalus, re The absorption cross 
sections for NoOz, HNO,, HNO, | CH, CO,, CH,O, and HO, are listed in 

Table B2. 

In computing the time dependence of atmospheric composition, the computation 
of the photodissociation rate coefficients using Eq. (54) can be a time consuming 
process, As the solar zenith angle changes, the column densities N. will change 
and, with them, the J,'s. It is desirable to have available tables of these coeffi- 
cients that are not dependent upon specific distributions of molecular oxygen and 
ozone. To kee, the tables perfectly general, therefore, the coefficients are com- 


puted and tabulated in Appendix B as functions of the column densities of Oo and O3. 


15, Hudson, R,D., and Mahle, S,H. (1972) Photodissociation rates of molecular 
oxygen in the mesosphere and lower thermosphere, J. Geophys. Res, 
77:2902-2914, a 

16, Ackerman, M. (1971) Ultraviolet solar radiation related to mesospheric pro- 
cesses, in Mesospheric Models and Related Experiments, G, Fiocco, FEd., 


ID. Reidel, Dordrecht, Netherlands. 
17. Heroux, L., and Swirbalus, R.A. (1976) Full-disc solar flares between 1230 
and 1940 A, J. Geophys. Res. 81:436-440, 
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OT ee 


the composition of the atmosphere is computed, the dissociation and ionization of 
atmospheric gases by incident solar radiation must be considered. In most calcu- 
lations, these photoprocesses are handled as chemical reactions, Consequently, 
the rate coefficients for dissociation and ionization of the various gases must be 
known, 

The photodissociation and photoionization rate coefficients for a gas particle at 
any altitude in the atmosphere is a product of the flux of the solar radiation arriv- 
ing at that altitude and its photodissociation, or photoionization cross section, 
Several excellent reviews of the absorption properties of atmospheric gases can be 
found in the literature. Among these are Hudson? and Johnston and Graham. AS 
Also, an excellent survey article by Turco!4 includes a discussion of photodisso- 
ciation processes for several gases, 

Beer's Law expresses the fractional transmission of radiation through a gas 


mixture 


= WA) . .-T(A) 
Ti.) = WA) = 51 
Io(A) OD 


where 1 (A) is the intensity of the incident radiation, I(\) is the intensity of the 
transmitted radiation, and T(A) is the optical depth of the absorbing layer, all at 
the wavelength A. The optical depth is given by 


TIA) = EG, QN, (52) 
1 


where a, (A) is the absorption cross section of the i'th type of particle and N, is the 
total number of particles of this type ina cm? column along the ray path to the sun. 
The rate cocfficient for the photoprocess for j-type particles at wavelength A 


is determined by 


, 


JA) = PG (A) L(A) TOA) (53) 


where 0;{A) is the cross section for photodissociation or photoionization of the gas. 


12, Hudson, R.D, (1971) Critical review of ultraviolet photoabsorption cross 
sections for molecules of astrophysical and aeronomic interest, Rev. 
Geophys, Space Phys. 9:;305-406, 

3, Johnston, H.5., and Graham, R. (1974) Photochemistry of NO, and HNO, 

n, J, Chem, 52:1415-1423, 


compounds, Can, J, 
14. Turco, R.P, (1975) Photodissociation rates in the atmosphere below 100km, 


Geophys. Surveys 2:153-192, 


—_ 
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Therefore, the average collision frequency for charge exchange is 


Vin = 1.73725 n ¥. (48) 


The average reduced charge exchange collision frequencies between like particles 
are given in Table 2, Again, it is assumed that the ion and neutral gas tempera- 


tures are equal, 


Table 2. The Average Reduced Charge Ex- 
change Collision Frequencies Between Like 
Particles 


10 1/2 


1.0x10 °°. H (2T) 


11 1/2 


1.6x10 °° 0 (2T) 


11 1/2 


1.7x10 °° N (2T) 


11 


3.0x1072 He (er)2/? 


1 1/2 


1.1x 1072? O, (27) 


11 1/2 


2.2x10- 


No (2T) 


The average ion-ion momentum transfer collision frequency is given by 


xo: 1 -1 
y.. = 1,3 n.{ ——— }sec (49) 
ij iia) 


assuming that the charge states of the colliding particles are both unity. 


The collision frequency between i and j-type particles is then simply 


mi 
v.. = ———_P,.. (50) 
(mj + mj) 4 


6. THE CHEMISTRY AND PHOTOCHEMISTRY 


The chemical reactions included in this model are listed in Appendix A. When 
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Table 1, Polarizabilities of 
Neutral Gases 


Let 
Qprae = (46) 


where Y is the chemical reaction rate constant and v is the velocity of the colliding 
particle, 


Since 


1/2 
Vv - (se) » then 
m 


1/2 
pn). ee (1) ae Soe eee (47) 
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5. COLLISION FREQUENCIES 


A good review of collision processes can be found in Banks and Kockarts, : 
All of the material that follows was taken from that source, 


The average neutral-neutral collision frequencies are given by 


1/2 
a Ss 2 [27kT =I 
vi 2n, 05; (2a) sec (43) 


where uw, the reduced mass, is given by 


Ves ar 
m. +m. 
1 J 
and 
O.. stig + g.) 
1 2 J 


where 0; and g, are the diameters of the colliding particles. 
The average collision frequency for momentum transfer between ions and 


neutral particles is 


1/2 
— ms -9 a -1 
ee 2.6x10 Cr) sec (44) 


where ny is the concentration of the neutral particles, @is the neutral gas 


atomic polarizability in units of ie. °* em, and wisinamu, Table 1 lists the po- 


larizability of several atmospheric gases, 9 


The average collision frequency for charge exchange between unlike particles 


is given by 


1/2 
Bin EP (#5) (any? Q (45) 


where Qr is the average charge transfer cross section for ion-neutral collisions, 


In Eq. (45), the ion and neutral temperatures are assumed to be equal. 
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The calculation of F K however, requires the upper boundary velocity at the 


end of the time step. At the upper boundary, oe NO’, O3, OH, H, HOp, H,0, 


H)O,, Ho, NO,, O5 (' Ag), and argon are all minor species. It is assumed, there- 


fore, that they are in diffusive equilibrium there. This means that their diffusion 
velocities are zero at the upper boundary, but they are still subject to the mean 
mass motion, 

The remaining species, however, do have diffusion velocities at the upper 
boundary that must be defined, It is impossible to calculate the upper boundary 
velocities from Eq. (18) because concentrations are required at K+ 1, Since this 
altitude is outside the volume, the concentrations are not known there. Therefore, 
to estimate the velocities at the upper boundary at the end of the time step, we in- 
vestigate the continuity equation, Eq. (3), at the upper boundary where vi = 0 


(since the upper boundary is outside of the turbulence region). 


On; ro] 
Bt = P, - nL; - dz (n;w;). (41) 


Setting the time derivative to zero and solving the resulting equation for 


dw, /dZ gives 


dw. w. dn. P. - n.L. 
esd og ok + id . (42) 
dZ n; dZ n; 


This ordinary differential equation can be solved numerically over the last height 
step to give values of Wi at the upper boundary. Since Eq. (42) is solved in every 
iteration, the final iteration will give the upper boundary value of Wi at the end of 
the time step. 

The new concentration profiles at t + At are therefore computed as follows: 
Starting at the lower boundary, the E and F arrays are computed in the increasing 
height direction up to the upper boundary. The new concentrations are then com- 
puted starting at the upper boundary and applying Eq. (27) in the decreasing altitude 
direction down to the lower boundary. This procedure is repeated over each time 
step until the solutions converge for that time step. Once a solution has been ac- 
cepted, the time step is doubled (provided it has not reached the maximum value set 
by the program), Using the concentrations and velocities computed at the end of 
the last time step, the equations are solved over the next time step, In this manner, 


the Solution advances, 
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2 Eq. (36) can be written as 
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.* 
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: ae 
S a 
t 

- . AZ 

= be. 7) [ mi K 
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‘ ee 
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r a 

Di. = De. (38) 
2 Eq. (28) can be rewritten at level K as 
3 xt 
; Ex - ~ ~h 
Lt [Bk- CKExK- 1] 
ce and 
ot is Weel ed ae 
hr Fy = = ; (39) 
oe K p44 
a K” “K“K- 1 
& 3 
oe It is obvious, then, that 
a a 
® 
" and 
= Dee. 
ee Fy = (40) 
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né+p* tly 
nf +t =o 9 (30) 
[1+ Lé +l aty 
E, and Fo at the lower boundary are simply 
E =O 
° 
and 
Fosnf*?, (31) 


Oo co) 


At the upper boundary, the time derivative of the velocities in Eq. (7) is set 


to zero: 


oe dn; (1+0;) dT 1 
mln, az * Taz tH |7~ Sry %- ¥. ia 


Solving Eq. (32) for dn, /dZ gives 


dn; rm 33 
— = ~-n. i > KT : Vij (w; = w,). ( ) 


Writing Eq. (33) at the upper boundary (K) in finite difference notation gives 


Pat eee cae oe 
K+ 1 K _  _£t1 r Nk i L+1 

= - = si "WwW: ‘ 34 

AZ x ne K a, ; Vij (ws Wik (34) 


Solving Eq. (34) for a ; gives 


£+1_ 4+1 m; 4Z £+1 
K 


Replacing bie, : in Eq, (25) written at level K with Eq. (35) gives 


m. AZ 
| ak [: . PAZ Fe Mn ois | ns ek} 


OEE aR TEs SRS gO A aR ES RS OR EES sate Sag eran chet BES sept Soil epee eR ade oe Ma tah alsa Big oe wh ey 
Be TOI Be ete eee FO Chetan g Sue. Sea Stedn a | SZ oe Ae 5) eek MESO rec Ne Maas Do gto piah ae Meet Ste Se ee ys een Pane tae Rtg Fe Mtge Deehds LR A Ie eet get in t ee 
sent BPs EP vate mT Maar Sage AE as aR Oe Tet AY de Ne 


Sag hcg ae Re Feige RS ee at 8 RS ETE Gi Bie MG NS gS NS ac Tigehe ual an ee Seah eR otek 9 «al, 
gk Bate Sgiest geen ee GR 6 gene OG) gel ae Re Pe Tg pt eee, ONS, ng igs eg oye ON eigen Oe OT nt as eee al Oe 
Fee Ras hile ee oe me tram oe, gee se Se ag Og SE we a gs 8 ae te Ee ee os as 


a From Eqs. (26) and (27), the following definitions are obtained: 


i. 
Ps 
bo At 
on E, = = 
[BR > Ch Bx - i| 

8 

= and 

= Lact 

m . goes ee ee Wd oa 

= et eae 

= [Br Cy Ey - 1] 

Boe Using Eqs. (27) and (28), the height profiles for the concentrations can be cal- 
culated provided né + land nf + 1 can be specified at the lower and upper bounda- 
ries respectively. The boundary conditions will be discussed in Section 4, | 

| 
It should be noted that only some of the terms on the right-hand side of the 
Loe equations for At, BY, c, and pt are specified at the beginning of the time step. 
2. The chemistry terms and the velocity sums are specified at the end of the time 
step. Because of this, the solutions require iteration over the time step. At the 
end of each iteration of the concentrations, the chemistry and the velocities from 
Eq. (18) are recomputed. New concentrations are then calculated using these up- 
= dated values. This procedure is repeated until, on successive iterations, the con- 
is centrations of all species change by less than 1 percent at all altitudes, 
& 

“y 1. BOUNDARY CONDITIONS 

eo At the lower boundary, the concentrations of No, Oy, H,0, Ho, and Ar are 

oe held fixed for all time at their initial values (n4 + 1 = constant), All of the other 

‘eo species in the transport mode prove to be chemically controlled at 50km, so their 

wl concentrations are computed from the simplified continuity equations, With the di- 

om vergence of the flux set to zero in Eq, (1), that equation becomes 

e dn, 

~ ees aS A (29) 

a If the derivative in Eq. (29) is replaced by simple differences, then, at the lower 

se boundary, we have 

é 15 


a in which 
oe 1 
at a3 22>. ee 
ee K [1+ At Ag Yj] 
28 
est S, = (b- 7, a2) 
oe > 
pecan — - - 
y Ss. a (1 I. Az). 
‘ 
oe It is possible, therefore, to write a tridiagonal matrix between the lower and 
ao 
aun upper boundaries for each species using Eq. (21). These tridiagonal matrices are | 
ae then solved for the vertical distribution of the concentration using the method of 
Por Richtmyer and Morton, !1 
“ Substituting the recurrence relation 
eee +1 +1 
° fe Bee eRe (23) 
ae into Eq. (21) gives 
’ at phtla ph stl _ at 4+1_ af i oaae 
ig A M+ 1 * By My Ck Fy - 1% Cera? Dee an) 
woe This can be written as 
iz ee ea! PGs £t+1 £ 
boo. : : ed 
a A Mea 1 FLB > Cy Ey - 1] my Big i ae (25) 
Solving Eq, (25) for ny ve gives 
be £ es 7 
os $44 Ax eA Det Si Pe 
eae k ~ L £ Ae d. 4 L . (26) 
- [Be Oe Pee i [Be ~ Oy By | 
e 
Oe oe Eq. (23) written at level k is 
£+1_ £+1 
ny = Eym, ,+ Fy. (27) 
11, Richtmyer, R.D., and Morton, K.W. (1967) Difference Methods for Initial- 


Value Problems, Interscience Publishers, New York, 
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both computed at level k, Also 


Q2 + ( v4)? 
1 ee 
Ai = aC ara tay na (17) 
Q2 sin? 1 + (F V:« 
i : ij) 
J 
The finite difference equations, Eqs. (13) and (14), can be solved for the velocities 


fae at the end of the time step: 


Wk 
a+ atl 
eae : 2 ¥T, At ime 7 MK 
We Re ae ane? weer We . ee L+1 ot [4% 
[1 + ata; Ey; m, AZ, ny 
+ At Zv.wit } (18) 
j 1) Ix 


This equation holds for both the ionic and neutral species using the above definitions 
of I and with d; = 1 for the neutral species. The total molecular and eddy diffu- 


sion flux at level k can then be written as 


rn ee es ere ee ae 1 ati2 KT, At 
ny Wy =O) 2 = 
[1 + 4ta; a m, AZ, 
J 
at+1 £+1 
net - ny (- r, 42)| 
ey Sea l> Bie 4 Ge et 
+ ny, At AW - M41 7 Dy 
j#i “Jk AZ, 
(l- TA 240} (19) 
where 
BP uo: 5 ct Og A e AI 
I T, 2° TES ae 


£+1 


A similar expression can be derived for @ _ ; at level k-1. 
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3. THE FINITE DIFFERENCE EQUATIONS 


If 1 denotes the grids on the time mesh and k the grids on the height mesh, 
then the continuity equations, Eq. (1), can be written in finite difference notation 
at the height level kas 


+14 Atl aati 

ny, Me pltl_ sti, sti an a5 

— _k. A Ly} - (12) 
At k k AZ, - 1 


where, for the sake of clarity, the subscript i has been dropped. Using the equa- 
tions of motion, Eqs. (7) and (9), the o,* 1 and ort ; will be replaced in Eq. (12) 
with their equivalent expressions, 

The equations of motion, Eq. (7), for the neutral species can be written in 


finite difference form at height level k as 


4+ 1 4+i1 £+1 
w - Ww kT n -n 
a See eee eS a ee et 
At m; + 
ijn, AZ, 
+ > Vi. witl (13) 
jFi J dk 
and Eq. (9) for the ionic species as 
4+ 1 & 4£+1 L+ 1 
Wye We Re Ae bP L+1 
mY h+ 1 uae z Wie ri LY; 
At mj ny AZ, j 
Z£+1 
+ : ; ; 14 
», "i Sk aa 


For the neutral species, 


(1 +a.) OT 1 


ao eh 
k T° 02 Hj; 


(15) 


and for the ionic species 


Sacra ee ee ec 
re + + . (16) 
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In computing the transmittance of ET'V radiation, O, Oy, and No are consid- 
ered as the primary absorbers, The photoionization rate coefficients are com- 
puted by using Eq. (54), where the sum is taken from 0,1 to 102.57nm, Radiation 
below 15nm is considered to be X-rays, and its ionization is partitioned among 
O, Oo; No, and N in accordance with Swider, 19 

The calculations also take into account the scattering of some radiation from 
the sunlit hemisphere into the nighttime hemisphere. These scattered radiations 
include Lyman @ at 121.57nm, Ly:nan # at 102.57nm, the HE I line at 58.4nm, 
and the HE II line at 30.4nm. For simplicity, the amount of radiation that scat- 
ters into the dark hemisphere is taken to be a fraction of the direct noontime radi- 
ation at the particular altitude under consideration. These fractions are: 

Lyman ao, 0,01; Lyman 8, 0.004; HE I, 0.0011; and HE I, 0.001, 
Finally, the contribution of cosmic rays to the photoionization rate coefficients 


is also included. This is simply“4 


a -17 
qa = 1320 n.. (55) 


6.3 Electron Impact Tonization 


Reactions 208 through 215 in the list of reactions in Appendix A are the ener- 
getic electron processes, The flux of energetic electrons appropriate to the cal- 
culations made here were supplied by Jasperse. 2° The ionization rate coefficients 


for the energetic electron reactions are computed from? 


1 = foléb(é)aé (56) 
§ 


where 0(€) is the electron impact ionization cross section and @(£) is the energetic 


electron flux, both at energy §, The ionization cross section is given by 


iss 


o(ty=f- S(,W) dW (57) 
I 


d 24, Nicolet, M., and Aikin, A.C, (1960) The formation of the D region of the 
ionosphere, J, Geophys, Res, 65:1469-1483, 


25, Jasperse, J,R, (1978) Private communication, 
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in which 
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where I is the ionization threshold for a particular ion state, q, = 6.51 x 10714 
ae em’, Cof, is a constant, W is the excitation threshold energy inev, and P, y, 
v, and Q are parameters obtained from the best fit of Eq. (58) to experimental and 
theoretical results, Table B22 lists the values of these parameters used in the 
model presented here, The values of the rate coefficients computed for the ener- 
getic electron reactions are given in Table B23, 


6.14 The Thermal Diffusion Factors 


The thermal diffusion factors are assumed to be zero except for the following 
Species: O(-0.27), On (+0. 08), H(-0. 39), H,(-0. 31), and He(-0.36), A discussion 
of the thermal diffusion factors can be found in Zimmerman and Keneshea, 26 from 


which source theSe values were taken, 


7. A MID-LATITUDE ATMOSPHERE SIMULATION 


The computer code presented here describes the distribution of 56 individual 
atmospheric species, Not all of them have the full coupled sets of equations in the 
system, however. Some gases in the atmosphere have such short chemical time 
constants that they will react chemically long before they can be moved out of a 
unit volume by any transport process. The vertical distribution of these species 
can be determined quite accuratcly by considering only the chemical reactions in 
which they are involved and ignoring any transport effects on them, 

The species for which the coupled sets of mass and momentum equations are 


solved are: 


O H N co AR HH 
Og Hy NZD) — COg HE* 
Ov.CAg) HO, Ne ot 
O3 120 NO O35 
OH H Qo NOg NO* 
N20 nt 


Those species whose distributions are determined neglecting transport are 


26, Zimmerman, S,P., and Keneshea, T.J,. (1976) The thermosphere in motion, 
J, Geophys, Res, 81:3187-3197, 
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oc'p), Nj; o*(?D), and the electrons, The height distributions of all 56 species 
are computed over the entire altitude range, which is from 50 to 400 km, 
- In addition, the vertical distributions of the following species are computed 
- over a reduced altitude range, The temporal variation of these species is also 
computed with the transport term removed from their continuity equations. The 
distributions of these species are computed between 50 and 120km because, above 
this altitude, they do not exist in significant quantities. These species are: 
Oo O; H,05.CO, HNO, 
OF NO, 03.H,0 NO, 
O35 H’.H,O NO".H,O | 
oO; H’. (H,0) NO’. (H,0) 
= 4 aia 2~'2 
a NO,  #H.(H,O) NO", (H,0) 
ca 2 2-°3 23 
a NO,  H’.(H,O) NO*.CO 
e 2-°4 2 
ie CO, “H;0".N, 
CO,  4H,0°.0H 
4 37° 
HO .C0, 
Neither negative cluster ions nor metallic ions are considered here, To ensure that 
the atmosphere remains electrically neutral, one of the negatively charged species 
is always computed from the imposed condition of balance of charge. Below 120km, 
x where there are several types of negative ions in addition to the electrons, the 
te most abundant negatively charged species is always computed from conservation of 
ee charge, Thatis, its concentration is computed as the difference between the sum 
cae of all the positive and the sum of all the other negative ions. Above 120km, where 
=. the only negatively charged particles are the electrons, their concentration is al- 
ae ways Set equal to the sum of all the positive ions. 
a 7.1 Initial Conditions 
e 


Throughout the development of the time-dependent solution of Eq. (1), the tem- 
peratures and the turbulent diffusion coefficients are held fixed at all altitudes in 
the code presented here, For the simulation developed in this report, both of these 


parameters were derived from data obtained during the ALADDIN I experiments of 


20 November 1970.2°:28 The temperatures shown in Figure 1 between 50 and 


125 km were taken from measurements of Theon and Horvath. a From 125 to 

160 km, the measurements of Golomb?" were used. Above 160km, the tempera- 

tures were taken from the 900 K exospheric temperature winter model of the U.S, 

Standard Atmosphere Supplements, 1966. A discussion of the turbulent diffusion 

coefficients obtained from the ALADDIN I experiments can be found in Keneshea et 

al. 8 The coefficients presented there represent the minimum turbulent diffusion 

coefficients determined from the observations on the chemical trails. Following 

arguments given there, the turbulent diffusion coefficients (Figure 2) used in these 
calculations are three times larger than those presented in that paper. : 

Keneshea et ar” found that stability considerations required that the height 
steps in the mesosphere not exceed 150m. Inthe thermosphere, however, stability 
could be maintained with somewhat larger height steps. Therefore, to reduce the 
amount of core storage required for saving the numerous arrays as well as to re- 
duce the amount of computer time per solution, a variable height step array was 
used, From the lower boundary, 50km, up to 90km, the height steps were fixed 
at 150m, From 90km up to the upper boundary, 400km, the height steps were 
computed by assuming that they were proportional to the pressure scale height of 
the atmosphere. This produced an altitude array with 692 grid points, with the 
maximum height step at the upper boundary of about 1.5 km, 

The initial concentration and average velocity profiles were taken from the 
final solution of an unpublished simulation using turbulent diffusion coefficients one 
third smaller than those shown in Figure 2, A discussion of a method for establish- 
ing initial concentrations on all the species is given in Keneshea et a8 for a situa- 
tion where no previous information is available on the vertical distribution of the 
various species, 

The solutions to Eqs, (1) and (18) are started at noon of solution day 1 with the 
above initial values. Because this is a real time-dependent calculation which takes 
into account the daily variation of the solar zenith angle, no truly steady state con- 
dition on the species concentrations is ever achieved. The calculations are con- 
tinued, therefore, until all species reproduce their concentrations diurnally at all 


27, Rosenberg, N.W., Golomb, D., Zimmerman, S.P., Vickery, W.K. nd 
Theony JS, (1973)-The ALADDIN Experiment — -Part 1, Dynamics, Space 
Research XI XUI, Akademie-Verlag, Berlin, 435- 


28, Philbrick, C.R., Narcisi, R.S., Good, R.E., face. H.S., Keneshea, 
T.J., Macleod, M.A,., Zimmerman, S.P., and Reinish, B.W. (1973) The 
AL ADDIN Experiment-—Part Il, Composition, Space } Research XII, 


Akademic-Verlag, Berlin, 441-448, 


20, Theon, J.S., and Horvath, J.J. (1972) ALADDIN neutral atmosphere measure - 
ments, Trans, Am, Geophys. Un, 533463, 
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altitudes to within 1 percent. This condition of diurnal reproducibility is generally 


arrived at after about 20 solution days, 


7.2 Thermospheric Results 


The results to be discussed here were taken from the final diurnal cycle of the 
calculations when the major and most of the minor species at all altitudes were 
varying by no more than 1 percent from values achieved 24 hours earlier, The 
vertical distribution of the major species Nos Oo, O, and argon are shown in Fig- 
ure 3, In this one-dimensional temperature-independent calculation, No, Oo, and 
argon show no diurnal change, having arrived at an equilibrium with the invariant 
temperature given in Figure 1, Also displayed in Figure 3 is the fixed profile of 
helium used in these calculatians, 

Because of the constant temperature profile used here, atomic oxygen also 
shows no diurnal variation in the thermosphere down to about 120km, Below this 
altitude, the chemistry dominates the production (dissociation of O,) and loss of 
this species. i Thus, we observe in Figures 3 and 4 large diurnal variations of 
atomic oxygen in the mesosphere, 

Presented in Figure 5 are the mesosphere and lower thermosphere noontime 
profiles of the remaining important neutral species, Carbon dioxide is invariant 
in the mesosphere and lower thermosphere up to about 140km. Above this altitude, 
the diurnal variation of the photodissociation process becomes important in the de- 
termination of the carbon dioxide concentrations, The variation of the CO,/Ar 
ratio calculated here is shown in Figure 6 and compares well to that determined 
from atmospheric measurements of Trinks and Fricke. an 

Isodensity contour plots of the minor ncutral species that exhibit a large diur- 
nal variation are shown in Figures 7 to 11. The numbers in all the isodensity plots 
presented here are the common logarithms of the concentrations, TheSe species, 
for the most part, are important in the determination of the electron density 
throughout the atmosphere, Above about 100 km, the charge distribution is domi- 
nated by NO; Os, om He’, and a, Thus, the resultant thermospheric charge 
balance follows the chain crcated primarily by the time-dependent distributions of 
the neutral species N, NO, and N(7D) as shown in Figures 7, 8, and 9, respectively. 
The other neutral thermospheric species in the code have little effect in determin-~ 
ing the charge distribution, 

The ion and electron isodensity contours determined from these neutrals are 


stioown in Figures 12 to 19, As expected, the dominant positive ion in the upper 


a0, Prinks, H., and Fricke, K.H, (1978) Carbon dioxide concentrations in the 
lower thermosphere, J, Geophys, Res, 83:3883-3886, 
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thermosphere is the Ot ion. The electron density distribution, therefore, follows 
the distribution of this ion very closely. Also displayed in Figure 19 is the altitude 
profile of f Fo: the space-time dependence of the peak of the electron density along 
with the measurements of this parameter by Evans, a4 The correspondence is 

quite good, clearly demonstrating the midday altitude and amplitude and the night- 
time transition to higher altitudes and lower amplitudes. No horizontal winds are 
included in this one-dimensional simulation; yet, we obtain an excellent represcnta- 
tion of the Fo layer maximum, 

Our rationale, then, is that the contribution of chemical production and loss 
coupled to the convergences and divergences created in our calculations are the 
main mechanisms that control the amplitude and altitude variation of the Fo layer 
maximum, The question that arises from this one-dimensional calculation, which 
apparently gives good agreement with some Fo layer maximum measurements, is: 
What is the effect of horizontal winds? The answer to this question may be in- 
ferred from measurements of Reinisch, 32 who determined the time dependence of 
fF for another period from ionosonde measurements, His results for fF 


o max 
are shown in Figure 20, We observe a smaller peak ionization amplitude than 


max 


Reinisch does, but, generally, the time dependence of his Foy layer maximum fol- 
lows our calculations quite well. Also observed are the large variations that pre- 
sumably arise from local winds that show multireversals during the measuring 
period. The point we wish to stress is that ionization calculated from internally 
determined neutral species distributions does result in an internally consistent Fo 
layer ionization behavior that compares well with some of the observed ionic beha- 
vior. An example of this is the recent work of Sethia et al, ae who have performed 
calculations similar to ours but include the effects of horizontal winds, They de- 
termine the distributions along magnetic field lines, while we calculate in the ver- 
tical direction with the magnetic field entering by way of the Larmor frequency and 
the magnetic dip angle. They input the neutral atmosphere, while we calculate it 
as an internally consistent integral part of the model, 

In Figure 21, we have superimposed upon our calculation of the total electron 
content that taken from Figure 3, model 4 (no winds) of Sethia et al. 33 We ob- 
serve that the shapes of the distributions are almost identical between the two 


one-dimensional models. The differences in magnitude can be attributed to the 


31, Evans, J,V. (1971) Observations of F-region vertical velocitics at Millstone 
Hill-I, Evidence for drifts due to expansion, contraction, and winds, 
Radio Sci, 6;609-626, 

32, Reinisch, B, (1983) Private communication, 

33, Sethia, G,C., Bailey, G.J., Moffett, R.J., and Hargreaves, J.K. (1983) 
The effects of neutral air winds on the electron content of the mid-latitude 
ionosphere and protonosphere in summer, Planet, Space Sci, 31:377-387, 
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different conditions under which the models were computed. We also compare 
(Figure 22) the total electron content to some measurements taken at Boulder by 


Sethia et ar as shown in their Figure 8 (July 2, no wind, and July 5, with wind). 
Again, the shape of our distribution compares well with the measured profiles. 
For those cases where the wind amplitude is not very large, the contribution of the 
Vx B mechanism is very much smaller than that of the chemistry, We make an- 
other comparison with the work of Poulter et al?4 shown in Figure 23, Here we 
examine the total electron content for measurements at Boulder and Lancaster 
(Figure 5 of Poulter et al), and again quite good agreement is seen between the 
shape of the distributions from the chemical-transport one-dimensional model and 
the measurements. The observed differences between the Boulder and the Lancas- 
ter curves result from the differences in the raypath geometries, Comparison is 
also quite good between our calculations and the {Fo measured by Kohl et a3° and 
discussed by Risbeth., Be 

We reiterate the point we wish to stress (see also Rush et a>"); A comprehen- 
sive, one-dimensional, time-dependent ion-neutral model demonstrates excellent 
agreement with measurements of electron densities. The inclusion of the effects of 
winds is necessary, however, to explain the deviation of the ions and/or electrons 
from their equilibrium distributions, which significantly vary over the diurnal 


cycle. 


7.3 Mesospheric Results 


The ionization in the mesosphere is much more complex than the ionization in 
the thermosphere. Ablating meteorites create long-lived atomic ions from the 
metallic species present in the atmosphere, These metallic ions are significant 
and are observed to be the dominant ions in sporadic E layers, Their chemistry, 
however, is not well known. For this reason and because large amounts of addi- 
tional computer memory and time would be required for their inclusion in the 
transport mode, they have been omitted from the calculations presented here. This 


omission might explain why our electron density is somewhat low in the 90 to 110km ‘ 


34, Poulter, E.M., Hargreaves, J.K., Bailey, G.J., and Moffet, R.J. (1981) 
Electron content modelling: The significance of protonospheric contents, 
Planet. Space Sci, 29:869-883, 


35, Kohl, H., King, J.W., and Eccles, D, (1969) An explanation of the magnetic 
declination effect in the ionospheric F-layer, J. Atmos, Terr. Phys. 
31:1011-1016, 


36, Rishbeth, H, (1972) Thermospheric winds and the F-region: A review, 
J. Atmos. Terr, Phys, 34:1-47, 
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37, Rush, C.M., Kempner, M.P., Anderson, D.N., Stewart, F.G., and 
Perry, J, (1983) Improving ionospheric maps using theoretically derived 
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Figure 3. The Time Invariant Height Profiles of No, Oo, 
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Figure 5, The Noontime Profiles of Hg, HgO, NO, N90, CO, 
and COg in the Mesosphere and Lower Thermosphere 
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Figure 6, The Ratio of Carbon Dioxide to Argon at Four Different 
Times During the Day. The curves show the diurnal variation of 
carbon dioxide in the lower thermosphere 
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Figure 8, The Iso-Density Contours for Nitric Oxide 
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Figure 33, Altitude Profiles of Electron Density From Three 
Rocket Experiments at Wallops Island by Mechtly and Smith?9 
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Figure 32. The Lso-Density Contours for NO’: (HyO), in the Mesosphere 
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Figure 29, The Iso-Density Contours for HO, in the Mesosphere 
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Figure 30, The Iso-Density Contours for Hi ,O. in the Mesosphere 
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Figure 27. The Iso-Density Contours for H,03 in the Mesosphere 
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Figure 23. The Total Electron Content. The continuous curve is 
for this model. The dashed and broken curves are the measure- 
ments of Pouiter et al. 
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region (Figure 19) compared to the work of Solomon et al. sa Below 90km, the 
cluster and negative ions change the effective electron density to a point where it 
disappears completely between 60 and 70km, 

The isodensity contours for the negative ions and the positive water cluster 
ions for the chemistry given in Appendix A are displayed in Figures 24 to 32. Iso- 
density plots for the remaining positive and negative mesospheric ions are not 
presented here because these ions never achieve a concentration greater than 
1 eee at any time during the diurnal cycle. Of course, the dominant negative 
ions in the mesosphere are apparently the cluster ions, Because their chemistry 
is still not fully understood, they have not been included in this work, 

Particular emphasis is given to the sunrise period and a comparison of our 
electron density profiles with sequential measurements made at Wallops Island by 
Mechtly and Smith. For clarity, the measured profiles are displayed in Fig- 
ure 33, and our calculated profiles are shown in Figure 34, In some instances, 
the comparison of theory with experiment is fairly good. Substantial differences, 
however, do exist. These are, in part, a result of the differences in the seasons 
(summer for the measurement and winter for the model). 

In explaining the D and lower E regions of the ionosphere, there is a question, 
still not completely answered, of the distribution of the positive and negative ions. 
To exemplify the puzzle, we cite two works, First, there is the report by Torkar 
and Freidrich of the measurements taken during the "Winter Anomaly" campaign 
in Spain, The comparison between our calculations and those measurements, 
shown in Figure 35, is surprisingly good, In contrast are the measurements of 
Arnold et al and Viggiano and Arnold, *7 These measurements, supposedly ex- 
plained by Viggiano et al, a show an extremely different distribution of the negative 
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Figure 19, The Iso-Density Contours for the Electrons, The 
crosses show the diurnal change in the altitude of the peak 
electron concentration in the model, The dashed curve indi- 
cates the measurement of this parameter by Evans 
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Figure 18, The Iso-Density Contours for N5. At night, the con- 
centration of Nb falls below lcm? at all altitudes 
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Appendix A 


The Chemical Reactions 


Table Al lists the complete set of chemical reactions used in this report. 
Reactions 1 through 127 have temperature invariant rate coefficients that are listed 
after each reaction along with the source from which the rate coefficient was taken, 
Reactions 123 through 184 are temperature dependent, The rate coefficients for 


these reactions are of the form 


| B 
k(z) = A (@)) exp aes, (Al) 


where the A's, B's, and C's are listed in that order after each reaction and T(z) is 
the temperature at altitude z, Reactions 185 through 207 are the photo-reactions, 
and are functions of the wavelength of the solar radiation and the local solar flux, 
Reactions 208 through 215 are the energetic electron processes, Some of the rate 
covtficients in the list of reactions have not been measured in the laboratory. 


Therefore, these rate coefficients are given as estimates, 
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Appendix B 


Photodissociation processes and Their 
Rate Coefficients 


Appendix B contains information on the photodissociation, photoionization, and 
energetic electron impact reactions, Table B1 lists the solar fluxes in the wave- 
length region 116, 2 to 730,0nm along with the absorption cross sections for Oo, 
03, HgO, HaO9, NO9, NO3, and N9QO, Table B2 contains the absorption cross 
sectioas for N90,, HNOg, HNOg3, CHy, COs, CH,O, and HOp. 

Table B3 lists the photodissociation processes along with their total photo- 
dissociation rate coefficient, The units of all the rate coefficients presented in 
this table and in all the tables that follow (Tables B4-B23) are sec}, The con- 
stants in Table B3 were computed from Eq, (52) using the entire unattenuated solar 
flux. The values listed here are indications of the maximum rate coefficient that 


each process can attain, 


BL. MOLECULAR ONYGEN 
Two distinct processes for O5 photodissociation are considered, 


0, hvy—~o+roctD), A <175 nm 


and 


QO, ’ hr —_ (@) + O, 175 < r < 242,2nm. 


. ov ony 4a ag ole eo ay Pe we Oe Beet Saeh edad, eo 3 : 
: Paes nes oo 2 . Bo Ta ok mar oe 3! a Sy sae As 

A : te é at a ot ele . = 4 7s ~ 
Sree ee eer OR, I, Ce ES a a a i a a) so Tete, ae 


16 


The cross sections for this molecule are taken from the compilation of Ackerman 

witti the exception of those waveleneths in the Schumann Range bands, The photo- 

dissociation rate coefficients for the first process above are listed in Table B4, 

and those for the second process are listed in Table B5, The column densities 

“iven in these and all the other tables that follow (Tables B6~ B23) are in units of 
5 


carn eoluim, 


B2. OZONE 
For ozone, two processes are also considered, 
-_ 1 Tes ‘ 
Q. ¢ hy om OCD) + OnC Ag), A 310 nm 
ary 


().} ht —™ OF O,, 310< A <730nm, 


’ 2? 


Oo absorption Cross sections for ozone are also taken from the compilation of 
‘ Woe Gs, ae Bock F 
“chopinan, The photodissociation rate coefficients for the process with the ex- 


cited products are listed in Table B6, and those for the process with the ground 


state products ace listed in Table B7, 


BS. WATER VAPOR 
Sobor vapor te phorodiesoctsiod primarily through the process 


nj}: hy we Off: TE, P2bl ee 2 188 anim, 


ab soorcifon oposs sections for fhis molecule are taken from Watanabe and 
katt, Pi photodbesociation vate coolfficients for this process are tabulated 
Pete V5 
ei he f 


and “olinory, ULI. (1993) Absorption coefficients of water 
Moria ci. vacnun ultraviolet, J. Opt. Soe, Am, 45:758-755, 


bt 


Bi. HYDROGEN PEROAIDE 


*% 


According to Urey ct ale and Holt ct al, Oy the photodissociation process for 


hydrogen peroxide is probably 
H,O, + hu-» OH + OH, A <300nm, 


The absorption cross sections above 2£5nm are from Urey et al, and those below 
225ntm are from Holt et al. The photodissociation rate coefficients for this molc- 


cule are presented in Table BO, 


BS. NIPROGEN DIONTDE 
Tin primary quantum yield for nitrogen dioxide is unity for the process 


NOY hr NO + O 


7 


for waveleneths from 398 to 244nin, The process 


NO, § hi» NO + ob) 


has unity primary quantum yield from 244 to 190nm, For these computations, 
bowever, all of the photodissociation of NO, is assumed to take place according to 


the later process, The absorption cross sections for this molecule are taken from 
_ Bt . BS ‘ ae : 
Nakasamaot al and from Hall and Blacet,. The photodissociation vate cocffi- 


cients for this process are listed in Vable BLO, 


Ro. THE NO: RADIOAL 
According to Jotimston and Graham, ~ the photolysis of NO., can proce as 


Pot dines, HlC,, Dawsey, T.f., and Rice, FLO, (1°29) The absorption spect our 
ainbdegetiposiion of hydrous peroxide by light, J. Am, Chom, Suc, 
Leet hs Lavon, 

| lott, ROB, Medbave, CORL, and Oldenbergs, O, (1948) Ultraviolet absorption 


Spectrum of avdvog. hn peroside, J. Chem, Phys, 16:225-229, 


Noasavama, Ul, Kitanvira, VIV., and Watanabe, K. (1959) Tonization pot orcas 


end Qosorplion ea ctiher nis ab tlivowen dioxide, J, Chem, Phys, aad dee: eke ea eal 
rol Ghadk, PeC., aed Bla ot, Bl) (1 92) Separation of the absorption sy. clipes 
Nog anu NO, tu the oan Jha Sno A, J. Chem, Phys, 20:174. 175 
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Table B3, 
Coefficients for Unattenuated Solar 


O2 + HV = 0 + 010 
O02 + HV = 0+90 
03 + HV = O10 + O210 


03 + HV = O + 02 
H20 + HV = OH + H 
H202 + HV = OH + OH 
NO2 + HV = NO + 010 
NO3 + HV = NO + 02 
N20 + HV = N2 + O10 


N205 + HV = NO2 + NOQ 
HNO2 + HV = OH + NO 
HNO3 + HV = OH + NO2 


CH4 + HV = CH3 + H 
cCO2 + HV = CO + 010 
CO2 + HV = CO + O 
CH20 + HV = CHO + H 
HO2 + HV » O02 + H 


Total Photodissociation 


1.67327£-06 
8.307SB8E-08 
9.72065E-03 
4.17995E-04 
7.39662£-06 
1.539387£-04 
1.51039E-02 
7.0A137E-02 
7.76036E-07 
6.92337E-04 
S.B47S6E-04 
1.16540£-04 
9. 08606E-06 
6.84043E-08 
8.24894E-09 
2.58902€-04 
1.13657€-03 


Rate 
Flux 


SEC-1 
SEC~1 
SEC-1 
SEC-1 
SEC~1 
SEC-1 
SEC-1 
SEC-1 
SEC-1 
SEC-1 
SEC~1 
SEC-1 
SEC-1 
SEC-1 
secm1 
SEC@1 
SEC=1 


Table B2. Absorption Cross Sections for N,O;, HNO,, HNO3, CH 


r 
COg, CHO, and HO, (Contd) 4 
WAVELENGTH BAND CROSS SECTICNS (CM2) 

ANGSTROMS N205 HNO2 HNOS CH4 co2 CH20 HO2 
3200.00-3250.00 6.00E-21 2.00E-20 14.00E~-22 0. 0. 2.20E-20 0. 
2250.00-23390.00 5.00E-21 2.00E-20 0. QO. QO. 3.40E-20 0. 
3300.00-2350.00 3.40E-21 3.00E-20 Oo. Q. QO. 1.S50E-20 0. 
3350.00-3600.00 3.00E-21 2.80E-20 09. 0. 0. 1.20E-20 0. 
3490. 09-2150.00 2.40E-21 3.10E-20 0, oO. 0. 2.00E-20 O, 
3450. u0--3500.00 2.00E-21 3.00E-20 0. QO. 0. 5.00E-21 0. 
3500.00-3550.00 1.50E-21 5.00E-20 0. 0. oO. 2.80E-21 0. 
355G.00-3600.00 1.20E-21 5.00E-20 9. QO. 0. 1.20E-21 0. 
3600.00-3650.00 9.00E-22 2.90E-20 0, 0. 0. oO. Oo. 
3650.00-3700.00 7.00E-22 8.00E-20 co. oO. 0. 0. oO. 
3700.00-3750.00 4.00E-22 4.00£-20 0, 0. 0. 0. 0. 
3750.00-3900.00 2.50E~-22 2.00E-20 oO. QO. oO. 0. oO. 
3390.00-3850.00 0. 4.00E-20 0, 0. 0. oO. 0. 
385C0.0C~3990.00 0. S.00E-20 0. 0. oO. 0. 0. 
3900.00-3950.00 O. QO. 0. 0. 0. Qo. oO. 
3950.00-4000.00 O. 0. oO. oO. 0. oO. oO. 
4000.00-4050.00 0. 0. oO. QO. oO. oO. oO. 
4059.90-4100.600 0. QO. oO. QO. oO. 0. oO. 
4190.00-4150.50 0. 0. oO. QO. oO. 0. QO. 
4150.06~4200.00 0. 0. O. QO. 0. oO. oO. 
4200.00-4250.00 0. oO. oO. 0. 0. 0. oO. 
4250.00-4300.00 0. Oo. QO. 0. 0. 0. QO. 
4200.00~4350.00 0. 0. oO. 0. QO. 0. oO. 
4350.00-4490.00 O. Oo. oO. QO. oO. 0. 0. 
4400.00-4450.00 0. Oo. oO. QO. oO. oO. oO. 
4450.00-4599.00 0. 0. 0. Q. 0. 0. oO. 
4500.00-4550.00 0. oO. QO. 0. 0. 0. QO. 
42°50, 00~6690.00 0. oO. Q. 0. 0. oO. o. 
4600.60-4650.00 0. oO. oO. 0. oO. 0. 0. 
4650.05-4700.00 0. 0. 0. QO. oO. oO. oO. 
4700.00-4750.00 O. QO. 0. Q. 0. QO. 0. 
4759,00-4800.00 0. oO. oO. QO. oO. oO. oO. 
AkUS.00-4850.00 C. 0. oO. oO. 0. oO. QO. 
$50.00-3900.09 0. 0. oO. oO. 0. 0. 0. 
490:).00-4959.00 0. oO. QO. QO. oO. 0. 0. 
4950.00-5000.00 0. 0. 0. oO. oO. oO. oO. 
9900.00-5050.00 0. QO. QO. QO. 0. 0. 0. 
5059.00-5100.00 0. oO. QO. QO. 0. oO. QO. 
5100.00-5150.90 O. QO. oO. Q. oO. 0. Qo. 
5150.00-5200.00 0. ° QO. QO. oO. oO. oO. 
%220.00-5259.00 0. 0. 0. Q. QO. 0. oO. 
+29,.90+-$390.00 0. ) QO. QO. QO. Oo. 0. 
~129.05-6359.00 0. 0. QO. QO. QO. 0. oO. 
~3509.90-5490.00 0. Oo. QO. 0. oO. oO. oO. 

tO. 09-£450.00 0. 0. oO. 0, oO. 0. 0. 
1y 80=-5007.90 0. Q. QO. oO. QO. oO. QO. 
xe dk) F0=$S850.00 0. (0) QO. QO. 0. 0. 0. 
me SE Die NO: P) Q. Oo. Oo. 0. oO. 

B23 ) Q. Q. Qo. QO. oO. 

Le OMG Oe fe) Qo. Oo. QO. 0. 0. 
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Table B2. 


WAVELENGTH BAND 
ANGSTROMS 


1792.60-1803.60 
1B803.60~-1816.40 
1816. 40-1830 .60 
1830. 60-1846.20 
1B46.20-1863.40 
1863. 40-1882.20 
1882.20-1902.40 
1902. 40-1924.00 
1924.00-1947.00 
1947.00-1971.80 
1971.80-1985.00 
1985. 00-2000.00 
2000. 00-2025.00 
2025. 00-2050.00 
2050. 00-2061 . 86 
2061 .86-2083.33 
2083. 33-2105.26 
2105. 26-2127 .66 
2127.66-2150.54 
2150.54-2173.91 
2173.91-2197.80 
2197.80-2222.22 
2222.22-2247.19 
2247. 19-2272.73 
2272. 73-2298 .85 
296.85-2325.58 
2325,58~-2352.94 
2352. 94-2380.95 
2380.95-2409.64 
2409. 64-2439.02 
2439.02-2469.14 
2469. 14-2500.00 
2500. 00-2531 .65 
2531.65-2564.10 
2564. 10-2597.40 
2597. 40-2631 .58 
2631 .58-2666.67 
2666.67-2702.70 
2702.70-2739.73 
2739.73-2777.78 
2777.78-2816.90 
2816.90-2857.14 
2857. 14-2896.55 
2898.55-2941.18 
2941.18-2985.07 
2985.07-3030.30 
3030.30-3076.92 
3076.92-3100.00 
3190.00-3150.00 
3°50.90-3209.00 


N205 


Qo. 
5.00E-18 
4.00E-18 
3.30E-18 
2.80E-18 
2.20E-18 
2.00E-18 
1,50E~18 
1.10E-18 
1.Q00E-18 
9.00E-19 
8.00E-19 
7.00€-19 
6.00E-19 
5.Q00E-19 
4.00€-19 
3.80E-19 
3.20E~-19 
3.00E-19 
2.40E&-19 
2.10E-19 
1.80E-19 
1.S50E-19 
1.30E-19 
1.00E-19 
8.00E-20 
6.00E-20 
4.00E-20 
3.30E-20 
2.80F-20 
2.10¢-20 
1.90E-20 
1.10€-20 
9.90E-21 


~-wonvooocooc00dda 


HNO2 


.QOE-21 
-5SO0E~21 
.QOE~21 
.00E-20 
-10E-20 


CROSS SECTIONS (CM2) 
CH4 


Absorption Cross Sections for N 
COo, CH,O, and HO, (Contd) 


90s, HNO,, HNO,, CH, 


HNO? co2 
Q. o. 1.30€-21 0. 
0. 0. 1.00E-21 0. 
0. 0. 6.00E-22 0. 
0. 0. 4.00E-22 0. 
0. QO. 3.00E-22 0. 
0. 0. 1.50E-22 0. 
0. 0. 8.00E-23 0. 
0. 0. 3.00E-23 0. 
0. oO. 2.00E-23 0. 
0. Qo. 1.00E-23 ©. 
8.00€-18 0. 8.00E-24 0. 
7.00E-18 0. 5.50€-24 0. 
5.00E-18 0. 4.00E-24 0. 
3.80E-18 0. 3.20E-24 0. 
3.00E-16 Q. 3.00E-24 0. 
2.40E-18 0. 3.00E-24 0. 
1.70E-18 0. 2.B80E-24 0. 
6.00€-19 0. 2.40E-24 0. 
4.60€-19 0. 2.00E&-24 O. 
3.40€-19 0. oO. QO. 
2.40€-19 0. Qo. 0. 
1.50€-19 0. oO. QO. 
1.10€-19 0. QO. 0. 
9.Q0E-20 0. oO. QO. 
8.00E-20 0. oO. oO. 
5.S50€-20 0. oO. oO. 
4.50E-20 0. 0. oO. 
3.80€-20 0. Oo. OQ. 
3.00€-20 0. 0. oO. 
2.80E€-20 0. oO. 3 
2.10E-20 0. oO. 3 
2.00€-20 O. oO. 4 
2.00E-20 0. QO. 6 
2.00E-20 0. QO. 7 
2.00E-20 0. oO. 8 
2.00€-20 0. oO. 9 
2.Q00E-~20 0. oO. 1 
2.00€-20 0. QO. 1 
1.80E-20 0. QO. 1 
1.60€-~-20 0. QO. 1 
1.20E-20 0. Qo. 1 
1.10€~20 0. QO. 2 
9.00E-21 0. 0. 2 
6.S0E-21 0. oO. 3 
5. 00E~21 0. oO. 3 
3.20E€-21 0. oO. 3 
2.30€-21 0. QO. 4 
1.80E~21 0. 0. 3 
7.Q0E~22 0. 0. 4 
3.90E~-22 0. oO. 3 


{ 


~) 


-SOE-2! 
-BOE-21 1.20£-18 
-OOE-21 1.C0E-18 
-OOE“21 5.00£-19 
-COE-21 3.B0E-19 
-OOE-21 2.60E-19 
-OOE-21 1.19€&-19 
-O0E-20 8.00E-20 
-20E-20 5.00E-20 
-SOE-20 2.20E-20 
-70E-20 
-90E-20 
-10E~20 
-50E-20 
-O0E~20 
-80E-20 
-O0E~20 
-OOE-20 
-50E-20 
-SOE-20 
-O0E-20 


CH20 HO2 


5.50€-18 
5.50E-18 
5.50E-18 
5.50E-18 
5.60£-18 
5.70E-18 
§.80E-18 
5.B0E-18 
5.80E-18 
5.80E-18 
5.B0E-18 
5.90E-18 
5.90E-18 
5.90E-18 
5.90€-18 
6.00E-16 
6.00E-18 
6.00E-18 
6.00E-18 
5.B0E-18 
5.50E-18 
5.00€-18 
4.S0E-18 
4,.00E-18 
3.70E-18 
3.00E~-18 
2.70E-18 
?.30E-18 
2.10E-18 
1.70E-18 


eogocoeo0o0u0090 


ve ee 


Table B2. Absorption Cross Sections for N,Oz, HNO,, HNO, CH, 
CO.,, CHO, and HO, 


WAVELENGTH BANO CROSS SECTIONS (CM2) 
ANGSTROMS N205 HNO2 HNOZ CH4 co2 CH20 O02 
] 
4 
: 121t5.67-1215.67 0O. QO. Q. 2.80E-17 7.30E-20 9.00E-18 0. 
1162.79-1169.59 0. 0. 0. 2.80E-17 1.10€-19 0. oO. 
1169.59-1176.47 0. 0. Q. 2.80E-17 8.80E-20 0. 0. 
1176. 47-1183.43 0. oO. 0. 2.80E-17 6.00E-20 O. oO. 
1183.43-1190.48 0. oO. oO. 2.80E-17 4.70£-20 0. 0. 
1190. 48-1197.60 0. 0. oO. 2.80E-17 3.80E-20 oO. 0. 
1197.60-1204.82 0. 0. Q. 2.80E-17 4.00E-20 0. Qo. 
1204.82-1212.12 0. QO. Qo. 2.80E-17 6.10£-20 0. 0. 
1212.12-1219.51 0. 0. 0. 2.80£-17 7.40€-20 0. Q. 
1219.51-1226.99 9. o. 0. 2.80E-17 1.00E-19 0. Oo. 
1226.99-1234.57 9. 0. 0. 2.80E-17 1.20E-19 0. 0. 
1234.57-1242.24 0. 0. QO. 2.80E-17 1.60E-19 0. 0. 
1242.28-1250.00 0. oO. QO. 2.80E-17 1.80E-19 0. 0. 
1250.00-1257.86 0. °. Q. 2.80E-17 2.10£-19 O. QO. 
1257.€a6-1265.82 0. oO. 0. 2-80E-17 2.80E-19 0. 0. 
1265.B82~1273.89 O. oO. 0. 2.80E-17 3.30E-19 O. QO. 
1273.89-1282.05 0. QO. 0. 2.80E-17 4,20E-19 0. Oo. 
4282.05-1290.32 0. 0. Q. 2.80E-17 5.90E-19 O. oO. 
1290.32-1298.70 0. oe Q. 2.80E-17 7.00E-19 0. 0. 
1298.70-1307.19 0. 0. Q. 2.47E-17 B8.00E-19 0. 0. 
1307.19+1315.79 0. QO. Qo. 2.24E-17 8.30E-19 0. Qo. 
1315.79-1324.50 0. 0. QO. 2-07E-17 8.50E-19 0. oO. 
1324.§0-1333.33 0. 0. oO. 1.57E-17 8.70€-19 0. Oo. 
1333.33-1342.28 0. oO. 0. 1.33€-17 8.70E-19 0. Oo. 
1342.28-1351.35 0. 0. Q. 1.20E-17 8.70E-19 O. QO. 
1351.35-1360.54 0. QO. QO. 1.00E-17 8.60E-19 0. QO. 
1340.54-1369.86 0. 0. QO. 6.40E-18 B.10E-19 0. oO. 
1369 86-1379.31 O. Q. QO. §.90E-18 7.20E-19 0. 0. 
1379.31-1388.89 0. QO. 0. 3.93E-18 6.40E-19 0. oO. 
13e8.89-1408.45 0. Q. QO. 2.24€-18 6.00E-19 0. 0. 
1408.45-1428.57 0. 0. 0. 6.72E-19 5.70E-19 1.00E-17 0. 
1429.57-1449.28 0. oO. oO. 6.72£-20 5.50E~-19 7.00€-18 0. 
1449,28-1470.59 0. 0. 0. 1.68€-20 5.70E-19 4.00E-18 0. 
1470.59-1492.54 0. oO. 0. 1.12E-20 5.70E-19 5.50E-186 O. 
1492.54-1515.15 0. 0. 0. 9.00E-21 S5.00E-19 1.20E-17 0. 
1515.15-1538.46 0. oO. 0. 5.60E-21 4,00E-19 1.80E-17 0. 
15§38.46-1562.50 0. QO. QO. 4.48E-21 3.00E-19 9.00E-18 0. 
1552.50-1587.30 0. 0. Qo. 2.80E-21 2.50E-19 B.00E-19 0. 
1587.30-1612.90 0. Oo. oO. 1.12E-21 1.60E-19 0. oO. 
1612.90- 1639.34 0. QO. °°. 8.00E-22 1.00E-19 O. 0. 
1639.34-1606.67 O. 0. oO. 0. 7.00E-20 0. 0. 
Yuo.67-16949.92 0. oO. 0. oO. 4.00E-20 1.50E&-18 O. 
1$94.92-1724.14 0. 0. 0. 0. 1.50€-20 1.20E-17 0. 
172d.14-1739.13 0. 0. 0. QO. 1.00E-20 6.00E-18 O. 
1699.93-1750.90 0. 0. QO. oO. B.O00E-21 1.80E-18 0. 
17539 39-1763.29 0. 0. 0. QO. 5.Q0E-21 0. 0. 
$°G5.20-1758.60 0. om 0. QO. 3.50E-21 0. om 
og HOet7TA 85 90. 0. QO. 0. 3.00€-21 0. oO. 
TAS ork Poet Os Oo. oO. Q. 2.2¢€-21 0. oO. 
Tol. oor tek nd 9. Qo. QO. QO. 1.5CE-21 0. 0. 


Table Bl, 


H,0,, NO,, NO,, and N,O (Contd) 


WAVELENGTH BAND 


ANGSTROMS 


5700.00-5750.00 
5750. 00-5800 .00 
5800.00-5850.00 
5850.00-5900.00 
5900,00-5950.00 
$950.00-6000.00 
6000. 00-6050.00 
6050.00-6100.00 
6100.00-6150.00 
6150.00-6200.00 
6200.00-6250.00 
6250.00-6300.00 
6300.00-6350.00 
6350.00~-6490.00 
6400.00~-6450.00 
6450.00-6500.00 
6500.00-6600.00 
6600.00-6700.00 
6700.00-6800.00 
6800.00-6900.00 
6900.00-7000.00 
7000.00-7100.00 
7100.00-7200.00 
7200.00-7300.00 


FLUX 


PH/CM2/S/A 


§.38E+13 
5.42E&+13 
$.42€+13 
§.42€+13 
5.44&+13 
5S.44E4+13 
S.42€+13 
5.40E+13 
5.40€+13 
S.40€&+13 
5.38€+13 
5.36€+13 
5.34€+13 
§.32€+13 
5.30€+13 
7.90E+13 
5.22€+13 
5.18&+13 
$.14E&4+13 
5.09F+13 
5.04E+13 
4.99E+13 
4.94E+13 
4.90E+13 


QO. 
oO. 
QO. 
oO. 
oO. 
oO. 
oO. 
QO. 
Q. 
oO. 
oO. 
0. 
QO. 


02 03 


4.75E-21 
4.55E-21 
4.35€-21 
4.42€-21 
4.61E€-21 
4.89E-21 
4.84€-21 
4.54E-21 
4.24E€-21 
3.90E-21 
3.60€-21 
3.43E-21 
3.17E-21 
2.74E-21 
2.61€-21 
2.40E-21 
2.07E-21 
1.72E~21 
1.37E-21 
V.V1E-21 
9.13E-22 
7.93E-22 
6 .40E-22 
S.146€-22 


eoooocooce 


soe 


Solar Fluxes and Absorption Cross Sections for O 


CROSS SECTIONS (Cm2) 


20 


H202 


NO2 NO3 


6.00E-19 
7.Q00E-19 
7.00E-19 
9.00E-19 
1.00E-18 
9.00E-19 
8.00E-19 
7.00E~-19 
5.00E-19 
4.00E-19 
2.00E-18 
1.80E-18 
6.00E-19 
4.00E-19 
2.20E-19 
2.00E-19 
3.00E-19 
2.00E-16 
2.0CE<-19 
QO. 
Q. 
QO. 
QO. 
oO. 


9+ Og. H 


29 | 


N20 


oO. 
QO. 
0. 
oO. 
QO. 
QO. 
QO. 
oO. 


om Table Bl, Solar Fluxes and Absorption Cross Sections for O 9» O3; Hy O, 
er WO NOo, NO., and N 90 (Contd) 
aera WAVELENGTH BAND FLUX CROSS SECTIONS (Cm2) 
ry) 4AGSTROMS P\\/CM2/S,a O02 03 H20 H202 NOo2 Nos n20 
cas 3200.00-3250.00 1.626+13 0. 1.50E-20 0. oO. 2.94E-19 0. QO. 
ne 3250.00-3300.00 1.94€&413 0, 7.78€-21 0. QO. 3.3B8E-19 0. 0. 
Mee 3300.00-3350.90 1.79E&+13 0, 3.72E-21 0. oO. 3.61E-19 0. oO. 
‘se 3359.00-3400.00 1.89E+13 0, 1.70E-21 0. 0. 4.17E-19 0, 0. 
eee 3400.00-3450.00 2.02E+13 0. 7.46E-22 0. 0. 4.30E-19 0. 0. 
——— 3450,00-5500.00 2.06E+13 0. 2.66€-22 0. QO. 4.48E-19 0. Oo. 
3500. 00-3550.00 2.06E+13 0, 1.09E-22 0. 0. 5.53E-19 0. 0. 
3550.00-36C00.00 2.08E+13 0, 5.49E-23 0. 0. 4.96E-19 0. 0. 
Os 3600.00-3650.00 2.36E+13 0. 0. QO. 0. 5.84E-19 0. 0. 
ae 2656.90-3700.00 2.46E+13 0. 0. oO. 0. 5.70E-19 0. Oo. 
} 3700.00-3750.00 2.48E+13 0. O. 0. oO. b S.74€-19 0. 0. 
' 3755.00-3A00.00 2.346413 0. 0. oO. 0. ' 6.52E-19 0. 0. 
bee. 2400.00~3850.00 2.22€413 0. 0. 0. 0. 6.45E-19 0. 0. 
3850.00-3900.00 2.18E+13 0. oO. oO. 0. 6.51E-19 0. oO. 
ie 3500.00-3950.00 2.38€&+13 0. 0. oO. 0. 6.08E-19 0. 0. 
39£0.00-4000.00 3.08E+13 0, oO. oO. 0. 6.68E-19 0. 0. 
, 4006.00-4950.00 3 BOE+13 0. oO. oO. 0. 6.52E-19 0. 0. 
4050.00-4:00.00 3.98E+13 0, 2.91€-23 0. Oo. 6.59E-19 0. oO. 
i 4100.90-4150.00 3.98E+13 oO. 3.14E-23 0. Oo. 6.03E-19 0. Oo. 
. 4150.00-4200.00 4.04&+13 0, 3.99E-23 0. oO. 6.24E-19 0. oO. 
my 4200.00-4250.00 4.02&+13 QO. 6.54E-23 0. 0. QO. 0. 0. 
elas 4250.00-4300.00 3.88E+13 0, 6.63E-23 0. 0. oO. oO. oO. 
4300.00-4350.00 3.96E+13 9, 8.66E-23 0. Oo. 0. 0. 0. 
4350.00-4400.00 4.50€+13 0. 1.25€-22 9. 0. 0. QO. oO. 
BeeO B05 1430-00 4.7B8E+13 0. 1.49E-22 0. oO. 0. Oo. Oo. 
50.0C~4500.00 4.96E+13 0, 1.71E-22 0. oO. Oo. oO. 0. 
ace 60-4550.00 4.98&+13 0. 2.12€-22 0. Oo. 0. 6.00E-20 0. 
4550.00-4600.00 4.96€+13 0. 3.57E-22 0. 0. Oo. 1.00E-19 0. 
4600.90-4650.00 5.00E+13 0. 3.68€-22 0. 0. Q. 8.00E-20 0. 
4659.00-4709.00 5.10E+13 0, 4.06E-22 0. Q. Q. 1.20E-19 0. 
479C.00-4750.00 5.22E+13 0. 4.89E-22 0. oO. QO. 1.60E-19 0. 
4750.00-4900.00 5.18E+13 0, 7.11€-22 0. Oo. Oo. 2.40€-19 0. 
\ 4800.00-4950.00 4.92E+13 0, 8.43E-22 0. Oo. 0. 2.190E-19 0. 
ane 4855 CO-4900.00 4.68E+13 0. 6.28E-22 0. 0. QO. 2.40€-19 0. 
r 4909.00-4950.00 5.06E4+13 0. 9.09€-22 0. oO. 0. 2.40€-19 0. 
pas 4950.00-5000.00 4.96E+13 0. 1.22E-21 0. 0. oO. 3.00€-19 0. 
Poets 5009.00-5050.00 4.98E+13 0. 1.62€-21 0, 0. 0. 3.00€-19 0. 
— 5020.00-5100.00 5.00E+13 0. 1.58E-21 0. 0. 0. 4.00E-19 0, 
bic 5100.00-5150.00 4.86&+13 0. 1.60E-21 0. oO. 0. 4.50E-19 0. 
. 5150.60-5200.00 4. 86€+13 0, 1.78€-21 0, Oo. oO. S.00E-19 0. 
: 52.30.00-5250.00 5.04E413 0, 2.07E-21 0. 0. oO. 5.00E-19 0. 
be 57 39.99-5390.00 5.16£+13 0. 2.556-21 0. 0. oO. 6.00E-19 0. 
aan 33v9.99-5350.00 5. 28E+13 0, 2.74€-21 0. Oo. 0. 6.00E-19 oO. 
Pgs £359.05-$400.00 5.34€413 0, 2.88E-21 0. Oo. 0. 5.00E-19 0. 
Laude b4t9.00 5.400413 0, 3.07E-21 0. 0. 0. 4.00E-19 0. 
oo SP a ue -S4y0U.00 $.30F+13 0, 3.17E721 0. oO. oO. 6.00E-19 0. 
Olle $550.00 5. 32F413 0, 3.36€-21 0. oO. oO. 6.50E-19 0. 
e : YoY Ld Slarlera a. 3.88E-21 0. oO. oO. 7.00E-19 0. 
feces Bog aw oN OGY. osko13 0, 4.31&-21 0 Qo. QO. 6.00E-19 Q. 
oe LuoS.G2-$700.00 5.34413 0, 4.67&-21 0 QO. QO. 6.00E-19 oO. 
(eM 
| 
e 
: 
6 74 
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Table Bl, Solar Fluxes and Absorption Cross Sections for Oo, O3, H,9O, HyO4 
NO, NO, and N,O (Contd) “ 


WAVELENGTH BAND FLUX CROSS SECTIONS (Cm2) 
ANGSTROMS PH/CM2/S/A 02 03 H20 H202 NO2 NO3 N20 
' 
f 1792,60-1803.60 %.01E+10 0. 7.77E-19 1.02E-18 0. O. 0. 0. 
= 1803.60-1816.40 1.27E+10 0. 7.60E-19 6.98E-19 0. 0. 0. oO. 
a 1816.40-1830.60 1.42E+10 0. 7.33€-19 3.64E-19 0. QO. 0. 0. 
a 1830.60-1846.20 1.29£&+10 0. 6.99&-19 1.80E-19 0. oO. QO. to) 
; 1826, 20-18663.40 1.31E&4+10 90. 6.56E-19 6.51E=-20 oO. oO. oO. oO. 
a 1663.40-1882.20 1.63E+10 0. 6.06E-19 0. Oo. QO. 0. Oo. 
¢ 1882.20-1902.40 1.92E+10 0. 5.45€-19 0. 9.79E-19 0. Oo. oO. 
‘ 1902,40-1924.00 2.29&+10 0. 4.86E-19 0. 8.10E-19 5.25E-19 0. 1.26E-19 
: 1924.00-1947.00 3.61€&+10 0. 4.36E-19 0. 7.03E-19 4.83E-19 0. 9.00E-20 
1947.00-1971.80 4.68&+10 0. 4.01E-19 0. 6.15&-19 4.00€-19 0. 8.Q0E-20 
1971.80-1985.00 6.21€+10 0. 3.42E-19 0. 5.64E-19 3.50E-19 0. 7.00E-20 
1485.00-2000.00 6.85E+10 0. 3.26€-19 0. 5.43€&-19 3.30E-19 0. 5 .SOE-20 
2000.00-2025.00 7.33E+10 0. 3.26€-19 0. 5.16E-19 3.50E-19 0. 4.00E--20 
1 2025.00-2050.00 8.74E+10 0. 3.28E~19 0. 4.79E-19 4.00E-19 0. 2.50E-20 
rey 20650.00-2061.86 1.10€4+11 1.05E-23 3.51E&-19 0. 4.66E-19 4.20E-19 0. 2.00E-20 
e 2061.86-2083.33 1.96€+11 1.00E-23 4.11&-19 0. 4.42E-19 4.30E-19 0. 1.70E-20 
< 2083.33-2105.26 3.33E+11 9.55E-24 4.84£-19 0. 4.22E&-19 4,40E-19 0. 1.10E-20 
2105. 26-2127.66 4.21£&+11 8.93E-24 6.26E-19 0. 3.92€-19 4.54E-19 0. 4.00E-21 
2127.66-2150.54 4.63E+1! &.28E-24 8.57E~19 0. 3.63E-19 4.58E~19 0. 2.00E-21 
2150.54-2173.91 5.73&+11 7.60E-24 1.17E-1E 0. 3.40E&-19 4.62E-19 0. 1.10£-21 
Pe 2173.91-2197.80 5$.S53E+11 6.92E-24 1.52£-186 0. 3.19€&-19 4.50E-19 0. 1,.00E-21 
a 2197.80-2222.22 7.08E+11 6.28E-24 14.97E-18 O. 2.98E-19 4.25E£-19 0. 4.00E-22 
2222.227-2297.19 7.23E+11 5.65E-24 2.55€-18 0. 2.67E~19 4.00E-19 0. 2.80E-22 
2247.19-2272.73 7.13E&+11 5.03E-24 3.24E&-18 0. 2.49€-19 3.50E-19 0. 2.00E-22 
. 2272.73-2293.85 8.65E&+11 &.40E&-24 4.00€-18 0. 2.36E€-19 3.00E-19 0. 1.20E-22 
4 2299.85-2325.58 &.98E+11 3.76E-24 4.83E-18 0. 2.16£-19 2.50E-19 O. 7.00E-23 
noe 2325.58-2352.94 8.22E+14 3.09E-24 5§.79E-18 0. 1.99E-19 2.00E-19 O. 5.00E-23 
: 2352.94-2380.95 7.89&+11 2.44E&-24 6.86E-18 0. 1,81E-19 %.18E-19 O. 3.00F~23 
2380.95-2409.64 8.09E+11 1.75E-24 7.97E-18 0. 1.66E-19 1.00E-19 0. 2.00E-23 
2499.64-2439.02 8.5t€&+1t 6.74E-25 9.00€&-18 O. 1.48E-19 7.00E-20 0. t.10&-23 
2439.02-2459.14 9.O07E+11 QO. 1.00€-17 O. 1.30E-19 5.00E-20 0. 9.00E-24 
‘) 2469.14-2500.00 9.33E+11 0. 1.07€-17 0. 1.22E-19 3.60E-20 0. 8.00E-24 
ti 2590.00-2531.65 9.54&+11 0 t.11E-17 0. 1.08E-19 2.50E-20 0. 7T.09E-24 
. 2531.65-2564.10 1.22E+12 0. 1.12€-17 0. 9.30E-20 1.50E-20 0. 6.0GE-24 
i 2554.10-2597.40 2.14€+12 0. 1.11E-17 0. 8.43E-20 1.70E-20 0. 5.0CE-24 
2597.40-2631.58 1.28&+12 0. 1.03E-17 0. 7.49£-20 1.90E-20 0. §.COrE-24 
.- 2631.58-2666.67 3.19€+t2 0. 9.43£-18 0. 6.39£-20 2.10E-20 O. 5.00E-24 
: 2666.67-2702.70 3.47€&+12 0. 8.23E-18 0. 5.64E&-20 2.90E-20 0. 5.00E-24 
2702.70-2739.73 3.13€+12 0. 6.81€-18 0. 5.00E-20 3.70E-20 0. 5.20E~-24 
2739.73-2777.78 3.13€+12 0. §.31E-18 0. 4.33&-20 4.20€-20 0. 5.906-24 
ps 2777.78-2816.90 3.53€+12 0. 3.99&-18 0. 3.62E-20 5.00E-20 O. 6.00&-74 
Ao 2B16.90-2857.14 4.226412 0. 2.84€-198 0. 3.27E-20 6.50E-20 O. 5.40£-24 
on 2957.16-2898.55 5.94E+12 0. 1.92E-16 0. 2.91£-20 8.00E-20 0. 5.20f-24 
2828.55-2941.18 9.15€+12 0. 1.14E-18 0. 2.53E-20 9.00£-20 0 5.00E-24 
27:341.18-2935.07 9.C9E+12 0. 6.60E-19 0. 2.15£-20 1.20E-19 0 5.10E-24 
3485.07-303C.30 6.536+12 0. 3.63¢-19 0. 1.75&-20 1.33E-19 O. 5.00E-24 
5°30.30-3076.92 1.09€+13 0. 1.97&-19 0. 0. 1.65E-19 0. S.206°74 
i @ Lcd 42-437 -3.00 2.578413 0. 1.257-19 0. Oo. 1.97E-19 0. B.40F -74 
poe dodo o.09 ti 2tiet3 0. 3.232-20 0. ©. 2.24E-19 0. 1.7GE-24 
5130.00-°3269.C0 1.39£+13 0. 2.31E-20 0. Qo. 2.61E-19 0. oO. 
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Table Bl, Solar Fluxes and Absorption Cross Sections for Oo, Oz, H,0, 


11,05, NO,, NO,, and NO = 
20°23 2 3 2 
WAVELENGTH BAND FLUX CROSS SECTIONS (CM2) 

ANGSTROMS PH/CM2/S/A 02 a3 H20 H202 NO2 NO? N20 
1215.677-1215.67 3.00&4+11 1.00E-20 2.32E-17 0. QO. QO. oO. QO. 
1162.79-1169.59 1.51&+07 2,00E-20 7.80£-18 QO. QO. Oo. oO. 0. 
1169.59-1176.47 3.87E&+07 1,25E-18 7.97E-18 QO. oO. 0. oO. 0. 
1176.4771183.43 1.61&+07 2.55E-19 8.66E-18 0. QO. oy 0. 0. 
1183.43~1190.48 1.76E&+07 3,00£-20 9.51E-18 9. 0. oO. Oo. 0. 
11S50.48-1197.60 2.55&+07 3.75E-19 1.25&-17 0. oO. QO. 0. 0. 
1197.60°1204.82 <«.87£+07 4.45E-18 1.84€-17 O. Q. 0. QO. 0. 
1204.82~-1212.12 1%.01E&+08 8,35E~-18 2.19€&-17 O. Qo. 0. QO. QO. 
1212.12~1219.51 3.08E+08 6.00E-19 2.30E-17 0. 0. QO. oO. QO. 
1219.5171226.99 4.91&+08 2.35E-19 2.26E-17 1.99E~-17 0. QO. 0. 0. 
1226.99~1234.57 1.80E&+08 4,.50E-19 2.06E-17 6.15E-18 0. 0. ; oO. oO. 
1234.57~1242.24 1.51&+08 3.35E-19 1.30€-17 5.16E-18 0. 0. 0. oO. 
1212.24~1250.00 9.69£+07 1.75&-17 8.91€-18 1.10E-17 0. QO. QO. 0. 
1259.00~1257.866 1.05€+08 B8.95E-19 7.24€-18 6.39E-18 0. Oo. -O. Oo. 
1257. 86-1265.82 1.19E&+08 4,30E-19 6.09£-18 7.13E-18 0. oO. 0. oO. 
1°55.02~1273.89 1.03&+08 1.10E-19 5.66€-18 7.84E-18 0. oO. 0. QO. 
1271.89~1282.05 7.96&+07 2,05E-19 S.87E-18 7.48E-18 0. QO. Qo. 0. 
1262.05~1290.22 5.71&+07 4.43E-19 6.47E-18 7.82E-18 0. QO. Qo. 0. 
1299.32~1238.70 7.16E6+07 5.55E-19 8.14E6-18 7.64E-18 0. QO. oO. 0. 
1298.70-1307.19 3.89&+08 4.20E-19 1.246&-17 6.71E-18 0. oO. oO. oO. 
1307.79~1315.79 2.06&+08 6.85E-19 1.52E&-17 6.75E-18 0. QO. oO. oO. 
1315.79-1324.50 6.68€+07 1,.45E-18 1.47£-17 5.67E-18 0. oO. QO. oO. 
1224.50~-1333.33 1.73€+08 2.25E-18 1.51&-17 4.93E-18 0. Oo. oO. QO. 
“123.339~1342.28 2.79E&+08 2.30E-18 1.51€-17 S.O0E-18 0. 0. Oo. QO. 
13462.28-1351.35 1.03€4+08 4,55E-18 1.65&-17 3.69E-18 0. QO. QO. oO. 
$351.35-1360.54 1.62€+08 7.23E-18 1.54E-17 3.36E-18 0. oO. 0. oO. 
1550.54~-1369.86 1.23E&+08 9.50E-18 1.35€-17 2.42&-18 0. oO. oO. oO. 
14€9.96-1379.31 1.37&+08 1,23€-17 1.05€-17 1.80&-18 0. QO. 0. 0. 
1379.31~-13388.89 1.20&+08 1.32E-17 7.97E-18 1.60E-18 0. QO. oO. 0. 

288.89-1408.45 2.93&+08 1.36E-17 7.17E&-18 1.19E-16 0. QO. 0. 0. 
1309.45-1428.57 2.10E+08 1,40&-17 6.28E-18 7.00E-19 O. 0. QO. 0. 
1424.57-1419.28 2.49£+08 1.48E-17 5.66E-18 S.O3E-19 0. 0. oO. QO. 
“sud. 28-1479.59 2.94€+08 1.41E&-17 5.23E-18 4.96E~-19 0. Qo. 8. 0. 
1379.59-1492.54 4.06€+08 1,29E&-17 4.47E-18 6.87E-19 0. 0. 0. 0. 
1492,.54-1515.15 4.55€&+08 1.15EF-17 3.69€&-18 9.4S5E-19 0. °. Oo. oO. 

“15.15-1538.46 5.96E€+08 9.91E-18 2.93E-18 1.38E-18 0. 0. Oo. 0. 
-539.46-1562.50 9.36E&+08 8.24E-18 2.19E&-18 1.84E-18 O. 0. oO. 0. 
1°52.50-15687.30 8.59E+09 6.58E-18 1.63E&-18 2.44E-18 0. QO. QO. 0. 
'587,30-°1612.90 8.05€+08 4.97£-18 1.20€-18 3.20E-18 QO. oO. oO. 0. 
$612.90-1639.34 1.29&+09 3.45E€-18 9.77E-19 3.91E-18 O. QO. QO. QO. 
12.°°9.274-1665.67 2.0€r+09 2.0dE-18 6.66E-19 4.69E-18 0. QO. 0. QO. 
'245.37-1694.92 2.864€+09 1,23€-18 9.1a5-1%5 4.57E-18 0. 0. QO. 0. 
174. 9241724.14 4.726+09 7.226719 8.17F-139 4.35E-18 0. on 0. 0. 
VIL dl t4-1799.13 4.998409 4,.56F-19 8.57E-19 3.57E-18 0. 0. oO. QO. 
12090 43-1750.00 59.49€+09 2.749E-19 38.40&-19 2.64&-18 0. 0. 0. oO. 
6°59. 0049752.20 6.252459 0 8.182719 2.649F&-18 0. QO. QO. Oo. 
1.03.9C-1753,60 T.LOCE+O9 O 8.04E-19 2.22E-18 0. 0. 0. QO. 

eR Qe AN eee BO y.0cE“*9 2.C4E-16 0. 0. QO. QO. 
Same Seg feet e680 J.45E-19 1.858-18 0. QO. oO. QO. 
Ce VP Tae Go p.orf+)3 0 7.892719 1.62E-'8 0. Qo. 0. Oo. 


CH,O + hy» CO +H, 


for wavelengths from 240 to 360nmn,. The absorption cross sections for this mol- 
Bl4 _ B15 
and Calvert et al. 


also vive the quantum efficiencies for these processes, Table B19 lists the tota! 


ecule are taken from McQuigg and Calvert These authors 
photodissociation rate coefficients for formaldehyde. If both of the above proces=.. 
are considered to occur for this molecule, then the values presented in Table B1! 
must be multiplied by the relative cfficiencies as given by McQuigg and CalvertBit 
and Calvert ct al819 for each of these reactions, 


B11. THE HYDROPERONYL RADICAL 
The primary photodissociation process for the hydroperoxyl radical is 
WO, +hy—OH+O A<456nm, 


The absorption cross sections used for HO9g are a composite of the measurements 
of Hochanadel et aor and Paukert and Johnston, pend The photodissociation rate 


coefficients for this radical are listed in Table B20, 


Bit. MeQuive, R.D., and Calvert, J.G, (1969) The photodecomposition of CHoG, 
CDeO, CHDO and CHoO-CDe20 mixtures at xenon flash lamp intensifics, 
J, Am, Chem, Soc, £1:1590-1599, 

Blo, Calves J. G5: Kerr, J,A., Demerjian, Soby.,-and MeQuigg, hp, 4172) 
Photolysis of formaldehyde as a hydrogen atom source in the lower atnios- 
phere, Sefence 175:751-752, 


Bis, Hochanadel, C.J., Ghormliey, A., and Ogren, P.J. (1972) Absorption srce 
trum and reaction kinetics of the HO, radical in the gas phase, J. ov 
Phys. 56:4426-4432, " 

Biv, Paukert, T,T., and Johnston, A.S, (1972) Spectra and kinetics of the lycra 


yD 


pevosyl free vadical in the gas phase, J. Chem, Phys, 50:2624-2630, 


Chom, 


~ ~~ re sae 


for wavelengths from116 to 164nm. Most of the CH, photodissociation occurs at 
ru Lyman da line 121,6nm, The values listed in Table B16 are the total photodis- 
-oclation rate coefficients for methane. If the above processes have equal effi- 
-isneies, then each process would have rate coefficients equal to half of the values 
.inowiin the table. The absorption cross sections for this molecule are taken 


frou. Watanabe et al. ae 


bil. CARBON DIOXIDE 


The photolysis of carbon dioxide can be proceed by two different processes: 


CO, thy CO+ OD) 116.2< A<166nm 


CO, + hy --CO+O 166<)A< 215nm, 


absorption cross sections for this molecule are taken from Inn et al, Bit 


iver, ite and Ogawa. pee The photodissociation rate coefficients for the pro- 
- vielding excited atomic oxygen are listed in Table B17, and those for the pro- 


. Yielding ground state atomic oxygen are listed in Table B18, 


ata. FORMALDEHY DE 


\ithough formaldchyde can be photodissociated in several ways, the primary 


“lytic processes are 


H,O + hy—» CHO +H 


Watanabe, K,, Zelikoff, M., andInn, F.C. Y,. (1953) Absorption Coefficient 
of Several Atmospheric Gases, Geophys. Res, Paper, > 2ec pt. 
APCRL-53-25, ae 

ii. fnn, FLC.Y., Watanabe, K., and Zelikoff, M. (1953) Absorption coefficients 
of gases in the vacuum ultraviolet, Part II, COy, J. Chem, Phys. 
2151648-1651, =~ 

Heimerl, J, (1°70) CO, absorption coefficient 1655-1825 A, J, Geophys. 


Res, 75:5574-5575, 


Ocawa, VM. (1071) Absorption cross sections of Og and COg continua in the 
Schumann and UV regions, J, Chem, Phys, 54:2550-2556, 
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fB9 


and Johnston and Graham. !° Its photodisso- 


cule are taken from Jones and Wul 


ciation rate coefficients are listed in Table B13. 


BY. NITROUS ACID 

The photodissociation of nitrous acid probably proceeds as 

HNO, + hy ~ OH + NO, 
The absorption cross sections for this molecule between 300 and 390nm are those 
given by Johnston and Graham, 13 Since O5 is not an efficient absorber of radiation 


in this wavelength interval, the photodissociation rate coefficients for HNO, are 
listed in Table B14 only as a function of the column density of ozone, 


BIO. NITRIC ACID 
The primary photolytic process for nitric acid is 
HNO. + hy — OH + NO, 
over the broad wavelength region from 192 to 325nm,. The values for the absorp- 


: ; : 1 dco 
tion cross sections for this gas are taken from Johnston and Graham, 2 The nitric 


acid photodissociation rate coefficients are listed in Table B15. 


Bil. METHANE 
The photodissociation processes for methane are 
rs hy — CH. +H 


arn 


Cll,» hype CH, ' Hy 


a 


me a 


Be, Jones, EVJ., and Wulf, O.R. (1937) The absorption coefficient of nitrogen 
pentoxide in the ultraviolet and the visible absorption spectrum of NO., 
J, Chem, Phys, 5:873-877, 


NO3 + hy-* NO + Os 


& 


for wavelengths longer than 578nm, For wavelengths bclow 578nm, it is not 

known whether the products are NOs + Oor NO+ Oo. In Table B11, all the photo- 
dissociation of this radical is assumed to produce NO + Op. The absorption cross 
sections for this molecule are taken from Johnston and Graham, i. Since rnolecular 
oxyven is not an efficient absorber for radiations above 45Unm, the photodissocia- 
‘ion rate coefficients for this molecule are given only as a function of the ozone 


coliumn density. 


7. NITROUS ONTDE 


The primary photolytic process for nitrous oxide according to Preston and 


2. 
wae ig 


NjO thy Ny + op), A<337nm., 


' a he : 1 ; 

‘low 210nm, oxygen atoms inthe S state may be produced. For the calculations 

presented in this report, Table B12 was computed with all of the photodissociation 
; 1 : : 

oy? NSO producing OCD). The absorption cross sections below 205nm are from 


“odioff et mee and above 205nm, from Bates and Hayes. Be 


BO. DINEPROGEN PENPONTDE 


Vho primary processes in the photolysis of NoO5 are not well understood, 


0 : ; 3 fee : 
Jones and Walt attvibute the pli otodissociation of this molecule to 


OO. PhyeN OO, + O-2NO, + O 


2 t 
roowevedenoths from 240 to s80am,. The absorption cross Sections fur tis tol 
e Saker 
Preston, Keo, and Barr, RRO (0°71) Primary processes in the photolysis 


of niteous oxide, J. Chem, Phys, vt os47-35-46, 


2, andinn, FPLC LY, (1955) Absorption co ficiom. 
oboeases inthe vactuin ulteaviolet. Part ll Nitrous oxide, J. Chom, Vhys, 


fotitoft, Mi, Watanabe, I 


) FAA eT ESET 
Mate <=, DLR rod Haws, POR. (les7) Atmospheric nitrous oxide, Vlanet, 
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Table B23. Rate Coefficients for Energetic Electron Impact Reactions 
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z(km) Solar Zenith Angle 
49, 78° 61. 50° 76. 22° 97. 16° 
95 0) 0 0 0 
100 3,54E-12 7, 83E-13 3.37E-15 0 
130 7.46E~-10 5. 79E-10 3.07E-10 0 
170 2.60E-09 2,25E-09 1.45E-09 0 
210 3.92E~-09 3. 69E -09 3. O3E-09 0 
250 +. 73E-09 4,.61E-09 4,22E-09 0 
300 3.27E-09 5.21E-09 5. 04E -09 0 
400 2,63E-09 2. 60E-09 2.52E-09 0 
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z(km) 
49. 73° 

95 0 
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130 4.11E-10 
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95 0 
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Appendix C 


The Computer Code 


The computer code used to develop the simulated atmosphere described in this 
report is presented in detail in this Appendix. 

To conserve computer memory, the code was written in OVERLAYs,. Once all 
the initialization is completed at the beginning of a run, there is no need to retain 
in memory all the code that accomplished this. 

The code is written in FORTRAN IV for a CDC 6600 mainframe computer, It 
also utilizes a large block of extended core storage for saving the concentration ar- 
rays, the average vclocity arrays, and the reaction rate coefficient arrays at the 
bevinninge of cach time step. If extended core is not available, these sections of the 
code have to be rewritten to use whatever disk or, as a last resort, tape storage is 
available, Using disk or tape storage significantly increases the residence time of 


the code on the computer. 


C1. THE WAIN OVEREAY (0.0) 


| 

Tic main overlay is relatively short. It defines all the blank COMMON vari- 

ebbeos shared by all the overlays. It also defines the parameters that describe the 
-invilation to be generated as well as the error criteria. It calls OVERLAYs (1,0) 


and (2,9), 
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C1. Blank COMMON 


There are 692 height steps between the lower and the upper boundaries in this 
model. The 424 height steps between the lower boundary and 120km represent the 
region over which the negative ions and the cluster ions are calculated. 

N is one of the dependent variable arrays. It contains the altitude distribution 
of the concentrations of the 26 species in the transport mode plus the arrays for 
oD), Nr, o* (2D), and the electrons. 

CONC is also a dependent variable array. It contains the height distributions 
for those species considered over the smaller altitude region. 

U is another dependent variable array. Here are stored the values at each al- 
titude of the average velocities of those species in the transport mode. 

NZERO and UZERO are working arrays that contain the values of N and U for 
the particular species whose set of equations is being solved at any given time. 

FORM and REMV contain the altitude arrays of the rate of chemical formation 
{P in Eq. (1)] and chemical removal [L in Eq. (1)] for the species whose equations 
are being solved. 

E and F contain the altitude arrays of E and F in Eq. (28) for the species under 
consideration. | 

Z is the altitude array in kilometers. 

DEDY is the array containing the fixed vertical turbulent diffusion coefficients. 


TEMP is the array containing the fixed temperature profile. 
DTEMD2Z is the array of the values 


1 dT (z) (C1) 


DELTAZ is the array oi height increments in centimeters, 

SUNSET is the array containing the sunset solar zenith angle at each altitude, 
GRAV is the sea level acceleration of gravity. 

COSD = cos (declination)*cos (latitude). 

SIND = sin (declination)*sin (latitude). Both COSD and SIND are used in com- 


puting the solar zenith angle, 


TIME is the independent variable. Time is tmevasured in seconds, with time = 
YQ as noontime of Solution day 1. 

CXI is the solar zenith angle in radians and is a function of the local time. 

RADIUS is the radius of the earth in kilometers. 

ISPECI is the total number of species included in the code. 

KSPECI is the number of species in the transport mode. 

NREAC is the number of temperature independent chemical reactions. 


NREAC2 is the number of temperature dependent chemical reactions. 
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K is the number of height steps over which the equations are to be integrated. 
KMIN1 = K - 1, 

ITURB is the index of the height of the turbopause in the Z array. 

CXINOON is the noontime solar zenith angle. 


SCXI is the local solar zenith angle in degrees. 


C1.2 Labeled COMMON 


Cl.2.1, COMMON/ESM/ 


OLDN, OLDCON, and OLDU are the extended core arrays corresponding to N, 
CONC, and U respectively. 

RATES is the extended core array containing the values of the height dependent 
reaction rate constants for the 88 temperature dependent reactions. These four 
Voriables are also defined in the LEVEL statement in all programs where they 


apply. 
C¥.2,2  COMMON/TIMES/ 


DA\5 is the variable that contains the number of solution days over which the 
‘quations will be solved. 

Dil. always contains the current value of the time step in seconds, 

{i MAX is the maximum value in seconds that the time step will be allowed to 
have, Th's value is 1800 seconds in code presented here. 

TSE? is the time in seconds after noontime of sunset at the lower boundary, 

TRISE is che time in seconds after noontime of sunrise at the upper boundary. 


(1.2.3 COMMON/ALTITUD/ 


vata ¥ ds the altitude of the lower boundary in kilometers, 

“TOP is the altitude of the u,-per boundarv in kilometers, 

YURDB is the altitude of the turbopause in kilometers, 

/7VINNS is the altitud: of the day-night transition height, During the day, it 


is sett Tian, and, during the night, it is set at 400 km, 
C12. COMMON/RESTART / 


LP RAAT is the number given to a particular solution file orn the interrupt file 


Tate $f) Vhose files are numbered sequentiallv from 1, 

TORE os the number of seconds that must be allowed the program to dump onto 
the dnteoreoot tile all the information needed to restart the code at the point when it 
sasocnt vaipted, The final solution of these sets of equations takes many hours of 
Ye ot! Computer time. Because of this, the code is normally run for a few 


fois tasualky 4), then stopped, and a dump is taken. Later, the code is restarted 


atepes port, continued for another period of time, and another dump is taken. 


This routine is repeated until the calculations reach a point where they repeat the 


solution over each diurnal cycle (usually about 20 days of problem time). 
TYME is the time in Seconds of the current computer run, The computer 


stops executing the code and dumps when 
TYME - (time running) < VOFF. 


TODAY is the date when the run was made and is printed out as XX/XX/XX, 
C1.2.5 COMMON/OUTPUTS/ 


TIMOUT is a variable that gives the code the frequency in ‘iours with which to 
send the solutions to the OUTPUT and the interrupt files. The solutions after 
every time step are not saved. Only those solutions occurring at TIMOUT hours 
and at interrupts are sent to OUTPUT and are preserved on Tape 4. 

OPRINT is the variable containing the irequency in hours at which solutions 
are to be written on Tape 4. Except on interrupts, TIMOUT is incremented by 
OPRINT every time a solution set is written to Tape 4. 

PRINT is a logical parameter under program control that determines whether 
writing to Tape 4 will take place after a solution is accepted (PRINT =. TRUE.) or 
will not (PRINT=. FALSE.). 


C1,2.6 COMMON/ERRORS/ 


ERRMAX is the criterion for acceptance of a Solution. When the relative 
change of all species at all altitudes over successive iterations is less than 
ERRMAX, the iteration is stopped and the solution for that time step accepted. 

EKSPECI is the criterion for stopping the iteration of individual species. 

Since the solutions converge faster for Some specics than for others, a good deal 
of computer time can be saved by no longer iterating the equations for those spe- 
cies whose relative change is less than FRSPFCI, Of course, ERSPECI must be 
much smaller than ERRMAX, 

ITMAN is the maximum number of iterations to be allowed for any given solu- 
tion. If any species have not converged after ITMAX iterations, the iterating is 
stopped and the last solutions computed are accepted, With ITMAX = 100, this sit- 
uation has never arisen, 

C1.2,7 COMMON/CONSTAN/ 

PT = 3, 14159265354, 

RADDEG = PI/180 is the conversion factor for degrees to radians. 

DEFGRAD = 180/PI is the conversion factor for radians to degrees. 

RADSEC = PI/43200 is the conversion factor for changing local time in seconds 


to local hour angle in radians. 
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BOLTZ = 1.38047 x 10718 ergs/K is Boltzmann's constant, 
SIN2I is the square of the sine of the magnetic dip angle. 
ATCON = 1, 660356 x 10724 is the mass of the hydrogen atom in grams. 


C1.2.8 COMMON/JAYS/ 


The labeled COMMON/JAYS/ contains the photodissociation rate coefficients 
for reactions 192 through 205, The locations reServed in this labeled COMMON 
are initialized from file Tape 3 by PROGRAM INITIAL, 


C1.2.9 COMMON/COLFRE/ 


The array D contains the invariant part of the collision frequencies Vj 5° 


NO} ed 


Lye 
Did, 1] = 5 (oe, + a,)° am (m, + “| 7 (C2) 
mymj 


Ci,2,109 COMMON/ATCONS/ 


The arrays in the labeled COMMON/ATCONS/ all pertain to the atomic and 
mole star constants. 

ATWT contains the atomic or molecular weights of the 26 species in the trans- 
port pode, 

MASS=ATWT*ATCON contains the mass in grams of the 26 species in the 
transport mode, 

Thi’ RV contains the thermal diffusion factors for the 26 species in the trans- 
mort imode, Where values are not known, the thermal diffusion factor is set to 


hich romoves any thermal diffusion effect for that species. 


“TET 
’ 


~YMBOL. contains the alphanumeric identification of each of the 56 species in 


Thee wneaecd, 

POLAR contains the polarizabilities of the species in the transport mode, For 
Speocheo so. cm polarizabilitics are not known, the value is set to zero to remove the 
polo. caotlita adbeet, 


Cl2.1) COMMON/RATECON/ 


tee arrays in labeled COMMON/RATECON/ all pertain to the rate coefficients 
stes ival reactions, 
srtass AK, BRK, and CK contain the A's, B's, and C's of Eq. (Al) for the 
Ste rebt copendent reaction rate coefficients. 
fie array DK contains the rate coefficients for the 215 chemical reactions at 
fora cre las altitude under consideration. 
areay CYB contains the alphanumeric cquivalence of the 215 chemical 


Poe tiomes themselves, 
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C1,2,12 COMMON/IONS/ 


The arrays in labeled COMMON /IONS/ all pertain to the ionization of the at- 
mospheric gases. 

WAVE contains the alphanumeric identification of the lines or wavelength bands 
in the EUV. 

FLUX contains the solar flux in photons /cm?/sec /A at the top of the atmos- 
phere at each of the above lines or in each of the wavelength bands. 

SIGMA contains the absorption and the ionization cross sections for Oo» No, 
and O, in that order. 

R contains the rate coefficients for the eight energetic electron reactions at 
the four zenith angles and eight heights listed in Table B22, 

The main overlay sets the values of the control parameters before calling 
OVERLAY (1,0). Upon return from OVERLAY (1,0), the initialization of the code 
for the runto be made is completed. The initial values of all of the dependent vari- 
ables N, CONC, and U are then saved in extended core storage. Finally, OVER- 


LAY (2,0) is called and retains control until the calculations are terminated. 


C2. OVERLAY (1.0) PROGR AVE ENEPIAL 


The purpose of OVERLAY (1,0) PROGRAM INITIAL is to initialize the code 
for the particular run to be made. In the dimension statement — 

S contains the collision cross sections in eine of each of the species in the 
transport mode, This information is required in the computation of the collision 
frequencies Vij: 

RAEN contains the ground level mixing ratios for the transport Species. 

SIG contains the photoici:zation cross sections for atomic nitrogen for the 
lincs and bands given inarray WAVE, These cross sections scaled by 1018 are 
listed in the DATA statement of this overlay. 

Following this DATA statement are the FUNCTION definitions. The first is 
ior the sea level value for the acceleration of gravity as a function of the latitude ®. 


[fis eiven ae 
A ’ bay a yo yop f ‘ ~.yeOOt : 2 = -6 . 2 ¢ 
GRAV(D) = 078.0556. 1 4 0,0052885 sin™ (6) - 5,9 x 10 sin (2@)]. (C3) 


The second is for the radius of the earth as a function of the colatitude and the sea 
level acceleration of gravity. It is given by} 


Cl. List, RiJd., Ed. (1963) Smithsonian Meteorological Tables, Sixth Revised 
iedition, Washington, D.C, 
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‘ 


altitude whose index is IALT. First, this subroutine computes the total concentra- 


tion Al as 


Then it moves the rate coefficients for the 215 reactions at the altitude w.ucr con- 
sideration into the K region, Following this, it ensures that the photodetachment 
and photodissociation of the negative ions are removed if the code is in the night - 


time. The remainder of the code computes the P and L for Species KIND and 


stores them in the appropriate FORM and REMV regions, respectively. The bulk 
or this subroutine as well as SUBROUTINE CHEMPR were not hand-coded for ob- 

viens reasons, To avoid the very tedious task of hand~-coding these routines every 
tit: the chemistry set or the species set is changed, a set of programs was writ- 
tor to perform the task of writing these codes. These programs can be found in 


. A od 
Rome shea, 


C3.) st BROUL TEESE RVUPEON 

SUBROUTINE RATECN sets the reaction rate coefficients for all the photo- 
veactions at the altitude whose index is IALT. The first three DATA statements in 
iiist subroutine contain, respectively, the solar fluxes (x io” photons), the ioniza- 
hay eross sections (x 10" cm”), and the effective absorption cross Sections 
(s ee cae for 0, (14g) as given by Huffman. 22 The next two DATA statements 
comtabip the altitudes (km) and the solar zenith angles corresponding to the input 
moos. or tie rate coefficients for the energetic electrons, Following this is the 
PATA statement specifying the rate coefficient for ionization by cosmic rays, The 
tiest oxccrtable statement of the program is the FUNCTION definition for the cor- 


‘rome nat altitude ZZ for the sea level acceleration of gravity. 


‘The eodo first computes all the rate coefficients for the current altitude and 
tore ihoan in the DN array. On the initial call to this subroutine, the fluxes and 
hee choss seetions viven in the DATA statements arc scaled to their proper orders 

Woot odiiude, and the vate coefficients for the energetic electron processes are 
reploood ds their natural logarithms, Next, the noontime values of H x F for Oo, 


ord QO are computed and saved for later use in computing the scattered radia- 


hho Tiere, THis the seale height of the species and F is the Chapman approxima-~ 

tron as fin computing the absorption, If the code is in the nighttime, then trans- 
rude ta statements near the end of the code since all photo processes are 

ae ee during the night, If the code is in the daytime, then the photo rate coeffi- 


ns are cietermined for the immediate daytime conditions, If the code is in the 


149 


fired sR SEM ecto AER 8 AS eae ihe Oe Wa igo BERS Ro tl hae raat Bias om oo ‘im! . ass ate i ad i =~ 


where, as before, £ and £+ 1 represent values at the beginning and at the end of 


the time step, respectively. The profile of o(4D) is computed once from Eq, (C21) 
each time that this subroutine is called, with Pe and oe used instead of p? + 1 and 
ioe i 2 Also computed in this subroutine are the concentrations of : (1) all the 
positive ions not in the transport mode; (2) all the negative species; (3) all the pos- 
itive cluster ions; and (4) HNO, and NO3 between 50 and 120km., 

Eq. (C21) requires the value of the chemistry terms F and L at the end of the 
timc step which in turn require the value of n at the end of the time step. Since 
this is the quantity being determined, Eq. (C21) will have to be iterated, Upon 
cach iteration, the chemistry terms are computed with the latest values of n, When 
the solution has converged, the value of n at £ + 1 will have been achieved and the 
chemistry will have been computed with this value, 

At the beginning of each iteration, the value of the concentration at the current 
altitucks is saved temporarily inthe array H, The new value is stored in the appro- 
priate cell of the F array. To ensure a more rapid convergence of the solution, 
the vali. of the concentration after the current iteration is replaced by the average 


vali: 


Thon the relative difference 


fon 
H 


is com: uted for cach species, If -he relative change for each species over the cur- 
root iteration is less than 20 percent, the solution is accepted for that height step 
and ihe iterations are started at the next higher altitude. If any species has a rela- 
five oo barmse creater than 20 percent, all the concentrations are iterated again. 
Phiie oro elute is vepeated until all species have relative changes at all altitudes of 
foes toap 24 pereent., The total number of iterations allowed at each altitude is 
ed, Phis rctimber has always heon more than sufficient to generate a solution, 

Soop th. concentrations of all the required species have been computed between 
this subroutine then computes the electron concentration between 


Peto and tie upper boundary from charge balance, 


(3.13 st BROL TINE CHEM 


~PRROTTPEINE CHEAL tompautes the quantities P and L for species KIND at the 
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The variable TIMOUT is used to keep account of the next hour in problem time 
when a solution will be sent to the output file. The code ensures that this time 
will occur exactly by resetting the time step if the current time exceeds this time. 
The RESET and PRINT switches are set to .TRUE., and the value of TIMOUT is 
incremented by OPRINT to update it for the next printout. When this routine is 
exited, DELT contains the time increment in seconds over which the next solution 
will be made, TIME contains the total elapsed probleni time in seconds, and CHI 


contains the solar zenith angle in radians at the end of the current time step. 


C3.b2 SUBROUTINE CHEMION 


SUBROUTINE CHEMION computes the concentrations of all species in the 
chemistry mode, This Subroutine also computes the largest negative species 
from charge balance. Below 120km, this could be a negative ion or the electrons. 
Above 120km, it is always the electrons. 

If the divergence of the flux in Eq. (1) is set to zero, then that equation can be 


written as 


dn. 
mode taki (C19) 
dt - Lol 


Pao convenient choices exist for solving these ordinary differential equations. 

They can be solved using the Kutta-Merson numerical integration scheme available 
in SUBROUTINE ODE, or they can be solved using simple finite differences. Since 
mich rmore computer time is required to solve the equations by the Kutta-Merson 
scheme and since it is generally a minor species whose concentration is being com- 
muted, the more economical finite differences are used to solve Eq. (C19). Writing 


Eq. (Cl') in finite difference notation at level k gives 


Ge Pog (C20) 
iS Leper th Get ty eral 
lt k k k 
wt . + 1 i 
Phe consentration ny at the end of the time step is then 
i £ 
ct] esa 
pot : (C21) 
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C3.8 SUBROUTINE ODE 


SUBROUTINE ODE solves a set of N ordinary differential equations using the 
Kutta- Merson numerical integration scheme. This program is explained in detail 
2 
elsewhere, The solution is taken over the altitude range DZ starting at altitude 


SZ. Yhe solutions are returned in array START. 


3.9 SUBROU PINE RHS 


SUBROUTINE RHS is called by SUBROUTINE ODE to supply to ODE the values 
of the right hand site of Eq. (42). The values are computed at altitude ZE using 
the intermediate values v computed in SUBROUTINE ODE. The results of the 
evaluation are returned in the array SLOPE. The intermediate values at altitude 
ZED are computed using the linear interpolation coefficients derived in SUBROUTINE 
VELOCITY, 


CBZ fost BROL EINE COLFRE 


SUBROUTINE COLFRE computes the sum of the collision frequencies (SNU) of 
species . sith all the other species at the altitude whose index is j. It also com- 
pLios the quantities = Ui w. required in Eq. (22) for species i. In addition, it sets 
the values oi omega (OMEG) in Eq. (17). OMEG is computed from Eq, (10) if i re- 


presents an ion, and is set equal to unity if i represents a neutral species. 


Cab SUBROUTINE FIMER 


Tixept for doubling of the time step after a successful integration and halving 
of dhe tins sfep on ancrror, all changes to the time and the time step are done in 
SUBROUTINE TIMER, This subroutine decides whether or not the next integration 
witl end ata print step, It also determines if the next time step will be in the day- 
tirae, techtiime, or in morning or evening twilight, The code advances normally in 
tite intii tais subroutine determines that the sun is about to rise at the upper 
Peoordars or is about to set at the lower boundary, When either of these situations 
is abort to nappen, this subroutine adjusts the time step so that the next integration 
will ond at the: beginning of morning or evening twilight. During these twilight pe- 
Wiods, the concentrations of some of the atmospheric species change rapidly. 
Those cnaneses fore the code to reduce the time step to successively smaller and 
sroalle co values before a solution can be accepted. To avoid the large amount of 
Con piter trike required to accompoish this, the code automatically reduces the 
they st cota do seconds with the onset of morning or evening twilight. The 
pew to Ad saiteh SPT controls this. The time increment is then allowed to double 


pn fe ROGRA ML INTER G after cach successful integration until it reaches TMAX, 
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If the ratio ns yy/n, > 100, the velocity is set to zero at that altitude. 


C3.7 SUBROUTINE VELOCITY 

SUBROUTINE VELOCITY sets up the parameters required for computing the 
upper boundary value of the velocity of species i. In order to use any numerical 
integration scheme to solve Eq. (42) for the velocity w at the upper boundary, some 


means of interpolating values of n, P, L, and dn/dz between the upper boundary 


and one step below must be provided. What is used here is a simple linear inter- 


polation: 
Nye = a2 +b 
(C17) 
OR ep eae oa 
solving these equations for a and b gives 
_"K  "K-1 
a = ———___—__. 
AZ. 
and 
Ne Zp 7 Nye Ape 8 
Ge! co Ne (C18) 
rae A ne 


inthe funetion defloitions at the beginning of SUBROUTINE VELOCITY, AF is 
the definition of a and BE is the definition of b. The subroutine then procedes to 
compute the aand b coefficients for the equations n(AN, BN), P(AP, BP), and L(AL, 
Glo. Phese cre stored in labeled COMMON/RHSIDE/ for use by SUBROUTINE RHS 
Ivete, Phis subroutine then cells SUBROUL.NE ODE to solve the ordinary differen- 


‘talequationsa  f.q, (42) for species i over the last height step. 


the total concentration 


the mean molecular weight in grams of the atmosphere 


(C15) 


eee 
m-=-=—, 
N 


and the mean mass velocity 


V2 (C16) 


The negative ions and the positive cluster ions, which are always very minor spe- 


cies, are excluded from the above sums. 


(3.6 SEBROUTINE EFORL 


For the species whose number is ICF, SUBROUTINE EFORU computes its E 
and F arrays of Eq. (28) if ICF is positive, or its average velocity array from 
tq, (18) if ICF is negative. 

The computation of the E and F arrays is fairly straightforward. At each alti- 
tude, the A(AJ), B(BJ), C(CJ), and D(DJ) values of Eq. (21) are computed. Then 
the i and I arrays are computed with the upper and lower boundary conditions im- 
posed, 

The specific velocities are also calculated in this subroutine since their evalua- 
tion requires much of the same coding as do the E and F arrays. Both computations 
oll SUBROUTINE VELOCITY to obtain the upper boundary velocity for those spe- 
vies that are not assumed to be in diffusive equilibrium there. To avoid unrealis- 
tically loree values of the velocities when the concentration is very small, the con- 
centr ition of each species n, and the ratio of n/n, are investigated at each alti- 
dey Of ny > 1.0 x 10 
dlues of VF < 1.0 x 10 


-3 J . 
em, then its velocity at altitude j is computed. For 


, the velocity is computed if the ratio 
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Performing the integration results in 


Mo (z ~ 2) 
TK | (2K 4 1) al 
R R 


where the subscript K refers to quantities at the upper boundary. After the concen- 
trations have been extrapolated to 500 km, the actual smoothing routine is called, 
where the profile is smoothed with a 7-point running mean. At this point, the code 
continues to loop back until all the required smoothing is accomplished. If no ion 
profile was smoothed, the subroutine is exited here. If at least one of the smoothed 


= Ny exp (C12) 


profiles was a positive ion profile, then the negative species must be readjusted to 
censure charge neutrality in the atmosphere. Below 120 km, the largest negative 


species igs recompuied from charge balance. 


n, = > n, - Dy; fie: (C13) 


“Sbove 120 km the electron concentration is recomputed from 


nar Dany. 


C2.bSEBROU TINE SMOOTH 


SUBROUTINE SMOOTH computes the smooth array B corresponding to the un- 
srpoothed array A. Both arrays contain N points. M, which must be odd, is the 
niraber of points in the running mean, In this smoothing routine, the beginning and 


The oond points care not changed, 


(asst BROEL LINE €O 
Ca adie cerre ponding te the index J, SUBROUTINE CO computes the 


density 


” a 4 at (C14) 
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assuming the concentration is in diffusive equilibrium above the upper boundary. 
That is, 


os + pg (C8) 
or 
dinkT) . _ amg (C7) 


which reduces to 


is 
= 
IF 


. __mg 
po dn? RE (C8) 


id 
nd 


It is also assumed that above 400 km 


dT . 
o— = 0. 


Therefore, Eq. (C8) becomes 


Tdn 2 | mg (C9) 
‘T 


ey 
jon 
N 
x 


An approximation to the acceleration of gravity at any altitude Z can be written 


g g .R™ 
ee cee ee) eee (C10) 


2 2 
(1+ %) (R +z) 


where ae is th » sea level value of the acceleration of gravity and R is the radius of 
the earth. Using Eq. (C10) and the assumption that the temperature is constant 


abave the upper boundary, the solution to hq. (C9) is 


n mgoR y d2-. 
“aq ny ee kT f (R + z)2 (C11) 
Zz 
fo) 
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have been computed for the species in the transport mode. It then calls SUBROU- 
TINE CHEMION to compute the concentration profiles at the end of the current time 
step of all the remaining species. After all the new concentrations have been com- 
puted, the code then calls SUBROUTINE EFORU to compute the average velocities 
for those species in the transport mode. 

Following this, a check is made to see if the solution for any species has con- 
verged at all altitudes. That is, is its largest relative change less than ERSPECI. 
If so, its IWHICH switch is set to zero and its ITER counter set to the current 
iteration number. The equations for this species will be skipped on the next itera- 
tion of the equation set. 

A check is then made to see if the largest relative change among all the species 
is less than ERRMAX. If so, the iteration loop is exited and the solution accepted 
for this time step. If not, transfer is made back to the beginning of the DO 130 loop 
for the next iteration. 

When the solution is accepted, the code then asks whether or not this solution 
should be sent to the output file. If printing has not been requested (PRINT= 
. FALSE, ), transfer is made to that section of the code where the time step is 
doubled if possible, and the new concentrations and velocities are saved in extended 
core storage. Transfer is then made to STATEMENT 1000, where a new solution 
is cenerated over the next time step. If the solution is to be sent to the output file, 
the code then calculates all the subsidiary parameters to be printed out along with 
the concentrations and the average velocities. The code at this point is fairly ex- 
plicit as to which quantities it is computing and sending to the output file. When all 
the information has been sent to the output file, the code then writes an interrupt 
file on Tape 4. This file contains all the information necessary to restart the 
code at this point at some later time. Since the development of a compiete time de- 
pendent solution takes many heurs of CDC 6600 CP time, it has to be accomplished 
in cuany short runs (about 4 hours each). ‘he code is restarted from Tape 4 each 


ti:ne. 


(3.5 SE BROL TINE SMOOTHY 


SMOOTHY controls the smoothing of those profiles that PROGRAM INTEG has 
determined should ne smoothed. If no profiles are to be smoothed (NUM=0), then 
the subroutine is exited immediately upon entry. Because of the nature of the 
smoothing process, it is neccessary to extrapolate the concentration profile above 
the upper boundary. The altitude above the upper boundary to which the extrapola- 
tion is dene should be far enough away so that the effects of the smoothing at the top 
oo tae apray are not felt at the upper boundary. In this code, the profiles are ex- 


trapolued to 500 km in l=km=height steps. The extrapolation is accomplished 
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Pe svlution, This may be a new solution over the next time step or another try using 
a smaller time step because an error was encountered. 

The code then saves the time, the solar zenith angle, the day-night transition 
height, and the day and night logical switches at the beginning of the time step. If 
an error is encountered later, these parameters will have to be reinitialized. The 
code then calls SUBROUTINE TIMER to compute the time and solar zenith angle at 
the end of the current time step. It then sets the JLOOK switch depending upon 
whether it is day or night. Following this, the IWHICH switches are all initialized 
to unity and the ITER locations to 999. Then, since the equations for helium will 
not be solved in this calculation, its IWHICH switch is set to zero. 

The statement DO 130 is the beginning of the iteration loop. The loop beginning 
with the next statement is taken over all the species in the transport mode. The 
equations will be solved over the altitude extent for each species in succession. On 


fae tirst pass through the loop, the quantities 


An Nan = Ne 3 
SS te 


uo e n AZ, 4 


(C5) 


ore co nputed and saved in the working array H. Next, the altitude dependent rate 
coofticients are computed using the current atmosphere and the values saved in ex- 
fond core storage. If the JILOOK switch is on.TRUE., then the photodissocia- 
lion rate cocfficients are sent to the output file. If any of the concentration pro- 
files aeed smoothing, this is done next and an acknowledgement is sent to the output 
Mile, To remain internally consistent, it is necessary now to recompute the velocity 
Praftife fos those species whose concentration profiles have been modified by 
smoothing, All the concentration- and velocities are then saved in extended core 
forage, 

‘The code then brings into memory the values of the concentrations and veloci- 
fies Tor the current species at the beginning of the time step. It then proceeds to 
Ceeabute oo new concentration profile in the DO 100 loop. It also saves the largest 
refative change for this species in the BiGEST region. The relative change is 


ePdMecdk cass 


Pho code fhen loops back to solve the continuity equation for the next species in 


morn. When the code exits at 105 CONTINUE, all the new concentration profiles 
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NSPEC is an array containing the order number from 1 to 56 of cach species 
whose concentration is to be smoothed before the iteration procedure is started. 
The variable NUM indicates how many profiles are to be smoothed. If NUM=0, no 
smoothing will be done. 

Two new labeled COMMONS are defined in PROGRAM INTEG. 


C321 Labeled COMMON/ LOGIC/ 


All the parameters listed in Labeled COMMON/LOGIC/ are logical variables. 
NIGHT is . TRUE. whenever the entire altitude region between the upper and lower 
boundaries is in darkness. Otherwise, it is .FALSE.. 

DAY is . TRUE. whenever the entire altitude region between the upper and 
lower boundaries is in daylight. Otherwise, it is .FALSE.. 

Both NIGHT and DAY variables are .FALSE,. when the sun is rising or setting 
and only part of the altitude region between the upper and lower boundaries is in 
darkness, 

JLOOK is set to. TRUE, to send to the output file, before the iteration begins, 
a table of the altitude distribution of the 05 concentrations, the 05 column densities, 
the 0% concentrations, the 04 column densities, and the rate coefficients for all the 
photodissociation reactions. If JLOOK is .FALSE., this table is not printed. 

RESET is set to. TRUE. in SUBROUTINE TIMER whenever the time step has 
heen vecomputed because the next solution will be sent to the output file and the 
interrupt file. 

tiND is set to. TRUE. whenever the computations for the current run are to be 


terminated, 


(3.2 Labeled COMMON /COEEL WM / 


The Labeled COMMON/COLUM/ contains information pertaining to the photo- 
dissociation rate coefficients. 

O2COL is used to transfer the 0, column density from SUBROUTINE DISSOC 
fo SEGROL TINE RATECN, : 

NICOL is used to transfer the 0., column density from SUBROUTINE DISSOC 
te st (Ot TINE RATECN, 

Din is the array where the photodissociation information is stored in SUBROU- 
TINE RATECN for printing in PROGRAM INTEG. To avoid excessively long print- 
outs, only every twelfth altitude value between the lower and the upper boundaries 
ts formed in these arrays. 

PROGRAAL INTEG starts by initializing some of its variables. STATEMENT 
ino CON PING Eo omarks the beginning of the working code. After every successful 


crunsite cessful solution, the code will loop back to this statement to begin the next 
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and (%, ‘This overlay uses two DIMENSION statements that contain the following 


variables: 

PHI is an array FQUIVALENCED to the array U(1, 21) and is used as tempo- 
rary storage of the total flux @ of each species during the printout phase of the 
code, The array U(1, 21) would normally contain the velocity profile of helium, 
but, since this species is held invariant in these calculations, these memory loca- 
tions are available for other uses, 

DED is an array FQUIVALENCED to the array E in COMMON. It is used as 
temporary storage of the turbulent diffusion velocity of each species successively 
during the printout phase of the code. 

DEDZ is an array EQUIVALENCED to the array NZERO and is used as tempor- 
tev storage of the vertical gradient of the turbulent velocity for each species during 
the printout phase of the code. 

VBAR is an array NQUIVALENCED to the array F and is used as temporary 
storage of the molecular diffusion velocity of each species during the printout phase 
of the code. The molecular diffusion velocity of a species is equal to its average 
velocity minus the mean mass velocity. 

Hi is an array used in various places in the code for temporary storage of 
arrays, 

OUT P contains four words each of which contains one hollerith character of 
the local time in hours. It is computed by SUBROUTINE HOUR. 

IWHICH is an array of integer switches, one for each species in the transport 
mode. If a switch is set to 1, the equations for that species will be iterated. If 
the switch is set to 0, the equations for that species will not be iterated. If IWHICH 
(21) is set to 0 before entering the iteration loop, the concentration of helium will 
remain fixed at its initial value throughout the calculations since its equations will 
not be iterated, 

IVER is an array that stores the number of iterations required by each of the 
scts of cauations before a final solution was obtained for a given time step. | 

RIGA is an array that contains the largest relative difference in each species | 


concentration profile. It contains the maximum value of 


ine appay ws 7 varies from the lower boundary to the upper boundary. This array 


is recomputed after each iteration. 
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given time. The first file, which is slightly different from all the others, must be 
prepared before this code is first run, It contains the initial conditions on all the 
dependent variables. (See the discussion on the initial conditions in Section 7.1,) 
Since this is an unformatted file, it must be written exactly as it is read here. 
First, it reads the date that the first file was written, the file number (1 in this 
ease), and the height in kilometers of the turbopause. It then skips over six words 
in the file that are not necessary to the current calculations. Then it reads the 
number of altitude steps in the complete arrays (K = 692), the number of altitude 
steps in the reduced arrays (K2 = 424), and the altitude array in kilometers. After 
this, it reads the height step array, the fixed temperature array, the array 

1/T dT/dZ, and the fixed turbulent diffusion coefficient array. Next, it reads the 
initial velocities and initial concentrations. 

In the following loop, the remainder of the interrupt files are skipped over and 
the tape is positioned at the beginning of a new file. At the completion of this loop, 
the locations containing the dependent variables will have been updated with the 
values cbtained after the last successful solution. These values are then used as 
initial conditions for the next time step. 

All files after file 1 on the interrupt tape are written in OVERLAY (2,0). Care 
must be taken not to write so many files on Tape 4 that it becomes full. Before 
this happens, a new Tape + must be made containing file 1 and the last file from 
the old Tape 4. 

After the dependent variables have been initialized, the code initializes the 
solar zenith angle to the current time and computes the noontime solar zenith 
anvlc, It then reads from input Tape 5 and, where required, conv rts to proper 
units the wavelengths, solar fluxes, and absorption and ionization cross sections 
10 te used in the photoionization calculations, Finally, from this same file, it 
reads the vate coefficient arrays for the energetic electron reactions, All that re- 
mains now is the calculation of the sunset (sunrise) solar zenith angles as a func- 
tion of altitude. Also calculated are the times of sunset at the lower boundary and 
of sunrise at the upper boundary. 

When all this initialization is complete, the concentration profiles correspond- 
ine to the initialized time are written to the output file. Finally, the current value 
of the variable TIMOUT is written to the output file in case the frequency of the 


ontpat of solutions is to be changed for the next run, 


(3, OVERDAY (2.0) PROGRAM INTEG 


OVERLAY (2,0) contains all the programs that actually solve Fqs. (1), (7), 
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RADUS (GRAV, @) 


2° GRAV 


2.27 x 10°? cos (2) - 2 x 10 7 cos (4) + 3.085426 x 10°° (C4) 

Next, the constants in labeled COMMON/CONSTAN/ are defined. After this, 
‘he constants for the particular solution to be generated are defined. 

CLAT is the latitude in degrees. 

Di'CL is the solar declination for the time of year when the calculations are to 
bo made, 

DIP is the magnetic dip angle appropriate to CLAT. 

Next, the acceleration of gravity, the radius of the earth, and the constants 
quired for later computing the solar zenith angles are computed. Information 
pertaining to the current calculations is then written to the output file (Tape 6). 

The next section of this OVERLAY sets the chemical symbols along with the 
aromie or molecular weight, the mixing ratio, the collision cross section, the po- 
larizability, and the thermal diffusion factor for the 26 transport species. Follow- 
ine this, the symbols for the species computed from chemistry only are written. 
Tne atomte or molecular weight of each species is then converted to mass in 
vranms, The ground level mixing ratio of No is computed to ensure that the sum of 
ail the mixing ratios is identically unity. The mixing ratios of the species are not 
howessars for the code described here, However, if other initial conditions are to 
In Consputed as discussed in Keneshea et al, : then the values of the mixing ratios 
ae roquleed, 

Tae following section of INITIAL obtains input information from various 

The list of chemical reactions with thei: rate coefficients must vu. pre- 
patron, separately and saved for us« here aS the tile called Tape 5, This informa- 
Pons coded according to the FORMAT defined in statements 770 or 771, depend - 
la cipor whether the reaction is not hcight dependent (770) or is height dependent 
(vd), As the reactions are read in, they are also sent to the OUTPUT file (Tape 6). 
Viooch voaction set is read in, the code reads in the photodissociation rate coef- 

-tatdes (Pables B4 through B20), All values in these tables are in log base 
lo, [ieve tables also are computed and saved Separately. Since these numbers 
nections only of the O. and O., column densities, they are universal. They 
rool be computed only once and can then be input to any calculation regardless of 
‘oruphic Location, season of the year, or time of day. 

Artop tine photodissociation rate coefficients are initialized, the invariant part 

© dlision frequencies y,, are computed and saved in the D array, 

(ie vest seetion of the code reads the interrupt file Tape 4, This tape is 


Stpactieod in files, each containing the solution for all the dependent variables ata 
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morning or evening twilight periods and the current altitude is below the day-night 
i. transition altitude, then the computation of the rate coefficients is skipped. Other- 
wise, they are computed. 
During the daytime, SUBROUTINE DISSOC is called for the computation of the 
a photodissociation rate coefficients. These are then stored in the proper cells of 
Ct the DK array. If the JLOOK switch is on .TRUE., these rate coefficients along 
with the concentrations and column densities of O5 and O03 are Stored in the DIS re- 
cion for later printing in the output file. Following this, the rate coefficients for 
the energetic electron reactions are computed, The values in the tables are inter- 
polated to the current solar zenith angle and altitude. Next, the column densities 
of O5, N,, and O at the current solar zenith angle are approximated. These col- 
wnaianeines are then used to compute the EUV attenuated flux at each wavelength. 
The rate coefficients for the photoionization processes are then computed. All 
@ livV radiation below 0.15nm is considered to be X radiation, which is handled dif- 
terentiv from the longer EUV radiation. The total x-ray ionization is computed and 
partitioned among the species according to Swider. 19 Finally, the total EUV rate 
coefficients are computed as well as the rate coefficients for the ionization of 
e nitric oxide by I.vman alpha and the ionization of ony. g). 
If the code is in the nighttime, the ahove computations are skipped and the 
code continues from this point. The contributions of the scattered radiations to 
the ionization rates are computed next, Finally, the contribution of cosmic rays is 


added into the rate coefficients and the subroutine is exited, 


(3.15 FE NECTION YING 


FUNCTION YINT performs interpolation where it is required. This function 
eascntiaily computes a cubic spline fit to the parameter under consideration, It 
thon returns the interpolated value in YINT corresponding to XII in the (X, Y) 


ena MS., 


(3.6 SUBROUTINE CHEAP VEAN 
SUBROUTINE: CHAPMAN computes an approximation to the Chapman function, 

which is roquired in computing the attenuation of solar energy in its transmission 

e through the atmosphere. The approximations computed are those given by Swider 
ous 

and Gardner, 
i} Be es 
. (7 oawider, We dre, and Gardner, M.k. (1967) On the Accuracy of Certain Ap- 
oy preximations for the Chapman Function, AFCRI.-67-0468, AD 658826, — 


C3.17 SUBROUTINE DISSOC 


SUBROUTINE DISSOC returns the value of the photodissociation rate coeffi- | 
cients appropriate to the current solar zenith angle and the altitude whose index is | 
IAL. '. The rate coefficients for all species except nitric oxide are determined as 
follows: The column densities of O05 and Og are first approximated for the current | 
O» and Oz concentration distributions in the code. A linear interpolation is then 
made in the appropriate photodissociation table to determine the current photodis - 
sociation rate coefficient for each of the photodissociation processes. Finally, 
using the data of Cieslik and Nicolet? as given in the DATA statement of this sub- 


routine, the photodissociation rate coefficient for nitric oxide is computed, 


C3.18 FE NCTION PRPEE2 

fUNCTION TRPLT2 interpolates in the two-dimensional photodissociation 
ratu coefficient tables. The two-dimensional array to be interpolated is fed to 
this function through the parameter ARRAY, The dimensions of ARRAY are (N,N). 
The tndices of elements in ARRAY are IO2 and 103. The increments required to 
pertoem the interpolation are DELO2 and DELOS. 


(3.19 FUNCTION PRPIEI 


FUNCTION TRPLT1 interpolates in the one-dimensional photodissociation 
rate cocfficient tables, The array to be interpolated is fed to the function through 
the pararieter ARRAY, The index in ARRAY is 103 and the interpolation increment 


is DRLOS, 


(3.20 56 BROUL TINE INDEN 


SUBROUTINI INDEN computes an index J] and an increment DEL corresponding 
to the column density COL. LIM is the lowest power of 10 of the column densitics 


of ©), and O., in the photodissociation tables under consideration. 


(3.21 St BROU TINE CHEVPR 


Pocte tine SUBROUTINE CHEMPR is called, it sends to the output file the 


complete list of chemical reactions for the altitude whose index is IZ, It also 


send= the rate coefficient for each reaction and its forward rate, This subroutine 
: F : 3 F 

iy also ccritten by the WRITER codes given in Keneshea.~ The output from this 

sabsoutine is used to determine the relative importance of the chemical reactions 


paoditPesent altitude regions, 


Ca, Cresiik, S., and Nicolet, M. (1°73) The aeronomic dissociation of nitric 
oxide, Planet. Space Set, 21:925-938. 
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oe C3.22 SUBROL TINE HOUR 

@ SUBROUTINE HOUR converts the value of TIME, which is in accumulated 
(se seconds from the noontime of the day on which the simulation was started, to the 
oe local time in hours (0000 to 2359), The digits of the hour are returned in hollerith 


in successive words of the OUTP array. 
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